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1.0  INTRODUCTION 


This  report  describes  the  systems  analysis  and  simulation  from 
the  start  of  the  XV- 5 A  Program  through  1  May  1963.  It  also  includes 
information  on  certain  related  topics  which  were  not  explicitly  a  part  of 
the  program.  This  task  was  performed  as  authorized  by  Contract  No. 
DA-44- 177-TC-715. 

The  analog  work  performed  included  the  six-degree-of-freedom 
hover  simulation;  the  longitudinal  three-degree- of-freedom  transition 
simulation;  the  longitudinal  three-degree-of-freedom  conversion  simu¬ 
lation,  the  six-degree-of-freedom  perturbation  conventional  flight  simu¬ 
lation,  and  the  simulation  of  gas  generator  control  for  wing-fan  thrust. 

Supporting  analyses  included  roll-yaw  coupling  and  structural 
feedback  in  the  pitch  mode. 

Further  work  accomplished  involved  the  stability  augmentation 
system  specification;  the  development  of  the  bridge  concept  for  roll  and 
yaw  louver  control;  support  of  various  hardware  tests;  the  generation  of 
the  specification  for  the  DeFlorez  point  light  source  visual  display; 
furnishing  consolation  services  during  the  DeFlorez  display  installation 
and  testing,  and  developing  the  yaw,  roll  and  pitch  direction  cosine 
relationships. 

This  report  is  divided  into  sections  corresponding  to  the  various 
subjects  mentioned  above. 
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2.0  SIX- DEGREE-OF- FREEDOM  HOVER  SIMULATION 


2. 1  Discussion 


The  hovering  simulation  utilized  a  six-degree  of  freedom  repres¬ 
entation  of  the  XV- 5 A  in  the  hovering  regime,  and  included  the  effects  of 
fan  ram  drag,  as  well  as  estimated  values  of  induced  lift  ahd  pitching 
moment. 

The  philosophy  of  using  only  rate  gyros  for  both  damping  and  ref¬ 
erence  was  carried  through  from  the  original  proposal.  The  original 
concept  is  shown  in  Figure  1.  A  separate  "Lagged- Rate  Integrator"  was 
used  as  a  short-term- attitude  reference. 


In  this  configuration,  the  time  constant  of  the  "approximate" 

integrator  was  chosen  as  10  seconds,  and  the  values  K  and  K^,  referring 

to  "rate"  and  "position"  gain,  were  independently  variable.  By  combining 

the  two  transfer  functions,  K  and  K  an  interesting  result  is  obtained. 

R  P 

1+rS, 


When  combining  transfer  functions  of  parallel  elements,  the  trans¬ 
fer  functions  of. the  original  items  are  simply  added,  thus  the  feedback 
terms  become: 


kr  +  kp 


1  +TS 


KR+KPtKRrS 
1  +TS 


(1) 


(Kp  +Kp)  (1  ♦  KRTS  )  (2) 


KRtKP 


1  +  tS 


Further,  letting  K  +  K  =  K,  and  lfetting  KR  =  R,  the  transfer 

KR+KP 


functions  becomes: 


(Xjj  +  Kp)  (1  +  -JL__)  =  K  a  +  RtS) 


“b 


+  K, 


1  +  rs 


(3) 


1  +rg 


Figure  1  Rate  Gyro  Concept  lot  Damping  and  Reference 


In  this  case,  the  system  gain  becomes  K,  and  R  is  defined  fes  the 
ratio  of  the  numerator  and  denominator  time  constants.  Since  R  Is 
always  <1,  this  transfer  function,  which  is  identical  to  the  original  para¬ 
llel  combination,  has  a  "lag-lead"  characteristic,  and  gives  the  effect 
for  short  periods  of  position  plus  rate,  when  fed  with  a  rate  gyro.  This 
network  is  easily  mechanized  with  two  resistors  and  one  capacitor.  Its 
use  is  discussed  in  another  section  of  this  report. 

The  hovering  simulation  was  utilized  to  give  a  preliminary  just¬ 
ification  for  the  lagged  rate  concept  as  applied  to  the  XV-  5A,  as  well  as 
to  the  selected  stability  augmentation  (SA)  system  authority  of  25  percent 
of  the  manual  input.  A  preliminary  investigation  of  handling  qualities  in 
the  hovering  regime  as  well  as  an  evaluation  of  hovering  thrust  control 
were  further  objectives  of  the  simulation. 

The  equations  of  motion  and  the  development  of  the  forces  used  in 
this  study  are  indicated  in  the  following  Section  2.3. 

All  evaluations  of  the  SA  system  and  hovering  flying  qualities  were 
made  using  a  fixed-base  cockpit  and  a  somewhat  marginal  display.  Past 
experience  indicates  that  results  obtained,  using  this  display,  have  been 
conservative,  since  various  vehicles  have  been  easier  to  fly  than  their 
simulation  indicated.  The  only  problem  inherent  in  a  marginal  display 
is  that  the  SA  system  can  be  "over-designed". 

All  of  the  hovering  analysis  and  simulation  used  an  assumed  trans¬ 
fer  function  for  the  electrohydraulic  servos.  This  transfer  function  was 
assumed  to  be  a  first  order  lag  with  a  1/10  second  time  constant, 
although  it  was  known  that  the  servo  would  be  at  least  third  order.  As 
far  as  aircraft  handling  qualities  are  concerned,  the  1/10  second  lag  is 
dominant,  and  the  higher  frequency  additional  second  order  response 
affects  only  the  servo  mode.  Later,  the  measured  servo  response  to 
electrical  inputs  was  determined  to  be  a  first  order  with  a  1/20  second 
time  constant  and  a  second  order  with  a  natural  frequency  of  20  cps.  As 
mentioned,  the  use  of  a  simple  first  order  lag  resulted  in  essentially 
equivalent  handling  qualities;  but  the  point  of  high-gain  system  instability 
was  not  accurate.  This  was  not  important  to  this  preliminary  investigation. 

Hovering  flights  were  made  to  evaluate  the  handling  qualities  of 
the  vehicle,  within  the  limitations  imposed  by  the  display. 

The  SA  system  includes  means  to  facilitate  maneuvering.  This 
involves  using  stick/pedal  motion  switches  which  modify  the  compensation 
network  whenever  the  control  inputs  exceed  a  certain  percent  of  full  scale 
deflection.  The  compensation  network  configuration  with  the  control  inputs 


5 


near  zero  is  called  the  "Holding"  configuration.  When  the  control  inputs 
are  displaced  beyond  a  given  amount  the  SA  system  is  in  the  "Maneuvering" 
configuration. 

In  the  holding  configuration,  the  compensation  network  provides 
the  approximate  attitude  stability  previously  described.  Both  the  ratio, 

R,  and  the  gain,  K^,  are  variable  in  the  primary  system.  In  the  maneu¬ 
vering  configuration,  simple  rate  damping  is  provided,  and  the  maneuver¬ 
ing  gain,  K  ,  is  variable  from  zero  up  to  3%  of  the  maximum  holding 
gain. 


As  a  result  of  these  investigations,  "nominal"  parameters  of 

about  K  =  20,  R  =  1  were  chosen  for  further  evaluation. 

n 

The  SA  system  as  presently  planned  has  authority  equal  to  25 
percent  of  the  manual  authority  in  the  various  control  axes.  It  was 
necessary  to  evaluate  the  effects  of  this  authority  limit  on  the  operation 
of  the  system.  The  authority  limit  has  two  effects.  These  are: 

1.  The  size  of  the  transient  generated  by  a  SA 
system  hard- over  failure  is  reduced  propor¬ 
tionally  to  the  reduction  in  authority. 

2.  The  ability  of  the  SA  system  to  recover  the 
vehicle  from  large  transients  is  also  reduced 
proportionally  to  authority. 

There  are  two  different  types  of  SA  system  hard-over  failures 
which  can  occur.  They  are: 

1.  The  amplifier  driving  the  SA  system  actuator 
Inputs  can  have  a  hard-over  output. 

2.  The  torque  motor  in  a  specific  actuator  can 
have  a  hard-over  failure. 

The  two  types  of  hard-over  failures  thus  have  different  results  on 
the  vehicle. 

An  amplifier  hard-over  in  any  axis  merely  results  in  a  25  percent 
maximum  displacement  of  that  axis  actuator(s).  This  difficulty  can  be 
overcome  by  either  switching  to  the  other  SA  system  channel  or  by  dis¬ 
placing  the  stick  to  restore  trim.  If  the  SA  system  is  not  switched  to  the 
other  channel,  the  axis  with  the  failure  will  have  no  SA  *  and  the  pilot 
control  task  becomes  more  difficult. 


A  torque-motor  failure  in  the  pitch  axis  results  in  a  complete  loss 
of  SA  in  the  pitch  axis,  and  the  stick  will  have  to  be  displaced  an  amount 
equal  to  the  authority  of  the  torque  motor.  Stick  inputs  will  still  operate 
the  pitch  buckets,  but  no  SA  is  available  on  primary  or  standby. 

A  single  torque-motor  hard-over  failure  in  a  wing  fan  louver 
actuator  results  in  a  1/8  full  stick  roll  mis-trim  and  a  simultaneous  1/16 
of  full  pedal  yaw  mis-trim.  The  mechanical  inputs  have  no  cross¬ 
coupling  induced  by  this  failure.  However,  the  SA  inputs  are  active  only 
in  3  of  the  4  torque  motors,  and  cross-coupling  results.  A  full  roll 
signal  would  result  in  3/4  of  the  normal  roll  authority  and  would  couple 
a  1/4  full  authority  yaw  signal  also.  The  opposite  would  be  the  same  for 
a  full  yaw  input. 

As  far  as  either  of  these  failure  modes  is  concerned,  the  25  per¬ 
cent  limitation  on  SA  torque  motor  authority  limits  the  effect  of  the 
transient  to  fairly  small  values.  It  is  easy  to  either  switch  to  the  stand¬ 
by  system  or  to  re-trim  the  vehicle  about  the  new  reference. 

The  second  result  of  the  SA  authority  limitation  involves  the 
magnitude  of  transients  which  the  SA  system  will  damp  out.  One  kind  of 
transient  which  is  easily  induced  is  derived  from  pushing  the  stick 
forward,  building  up  a  translational  velocity,  and  then  releasing  the 
stick.  This  is  not  an  expected  transient,  because  it  is  not  likely  the 
pilot  would  ever  do  this.  It  does,  however,  indicate  the  efficacy  of  the 
25  percent  authority  limitation. 

It  was  determined  that  the  SA  system  configuration  used  with  25 
percent  authority  would  take  care  of  a  transient  of  this  nature  up  to  +45 

feet  /second  forward  velocity. 


The  effect  of  possible  de- stabilizing  ground  effects  was  also 
considered.  The  worst  ground  effect  would  be  to  cause  a  first-order 
divergence  in  roll.  If  this  divergence  built  up  to  double-amplitude  too 
rapidly,  it  would  be  beyond  pilot  capability  to  control  it. 

The  SA  system  cannot  eliminate  a  ground  effect  of  this  type,  but 
it  can  raise  the  time  constant  to  an  arbitrarily  high  value,  and  thus 
reduce  the  divergence  to  a  controllable  value. 

The  reason  the  chosen  SA  compensation  network  cannot  remove 
a  pure  divergence  is  due  to  the  fact  that  it  is  a  rate  system,  and  this  leads 
to  a  zero  at  the  origin  on  the  complex  plane.  A  pure  divergence  is  repre- 


sented  by  a  pole  on  the  positive  real  axis.  From  root  locus  analysis,  as 
the  system  gain  is  increased,  the  roots  of  the  closed  loop  move  from  the 
poles  of  the  open  loop  transfer  function  toward  the  zeros  of  the  open  loop 
transfer  function.  Since  the  divergent  real  root  moves  along  the  real  axis 
toward  the  zero  from  the  rate  gyro  at  the  origin,  it  is  seen  that  the  diver¬ 
gent  root  can  never  migrate  to  the  stable  portion  of  the  complex  plane.  It 
can  only  be  moved  to  an  arbitrarily  small  distance  from  the  origin,  along 
the  positive  real  axis.  The  time  to  double-amplitude  of  the  response  of  a 
positive  real  root  equals  .  693/  a ,  where  a  is  the  location  of  the  root  on 
the  positive  X  axis.  Thus,  the  smaller  the  value  of  a ,  the  longer  the 
time  to  double-amplitude  becomes. 


t'Y’j’ 


2.2 


List  of  Symbols 

m  Vehicle  mass,  slugs 

m  Flow  rate,  slugs/sec 

I  Moment  of  inertia  about  the  body  X-axis,  slug  ft. 

x 

I  Moment  of  inertia  about  the  body  Y-axis,  slug  ft. 

y 

I  Moment  of  inertia  about  the  body  Z-axis,  slug  ft. 

z 

2 

I  Product  of  inertia,  slug  ft. 

xz 

H  Component  of  engine  and  fan  angular  momentum 

about  the  body  X-axis,  slug-ft.  2/sec. 


2 

2 

2 


H 

z 


S, 


l 

z 


Component  of  engine  and  fan  angular  momentum 
about  the  body  Z-axis,  slug  ft.  2/sec. 

/  2 

Acceleration  of  gravity,  ft.  /sec. 

Direction  cosine  between  the  inertial  Z  and  body 
X-axes. 


m 

z 


Direction  cosine  between  the  inertial  Z  and  body 
Y-axes. 


n 

z 


Direction  cosine  between  the  inertial  Z  and  body 
Z-axes. 


x  Distance  along  body  X-axis  (positive  forward  of 

the  center  of  gravity),  ft. 

y  Distance  along  body  Y-axis  (positive  out  the  right 

wing),  ft. 

z  Distance  along  body  Z-axis  (positive  down),  ft. 

h  Altitude  (measured  along  inertial  Z-axis),  ft. 

mf  Subscript,  denotes  main  fans. 
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Subscript,  denotes  left  main  fan. 

Subscript,  denotes  right  main  fan. 

Subscript,  denotes  nose  fan. 

Subscript,  denotes  commanded  value. 

Thrust,  lbs. 

Main  fan  thrust  for  rpm,  zero  vector  and  stagger 
angle  and  zero  speed,  lb 

Component  of  main  fan  thrust  along  the  body 
X-axis,  lbs. 

Component  of  main  fan  thrust  along  the  body 
Z-axis,  lbs. 

Nose  fan  thrust  (only  along  the  body  Z-axis),  lbs. 

Rolling  moment  (positive  for  right  wing  down), 
ft.  -lbs. 

Pitching  moment  (positive  for  nose  up) , 
ft.  -lbs. 


Yawing  moment  (positive  for  right  wing  back), 
ft.  -lbs. 


Force  along  body  X-axis  (positive  aft),  lbs. 

Force  along  body  Y-axis  (positive  out  right  wing) 
lbs. 

Force  along  body  Z-axis  (positive  down),  lbs. 

Velocity  along  body  X-axis,  ft.  /sec. 

Velocity  along  body  Y-axis,  ft.  /sec. 

Velocity  along  body  Z-axis,  ft.  /sec. 

Angular  velocity  about  the  body  X-axis,  rad.  / 
sec. 

Total  velocity,  fps 


Angular  velocity  about  the  body  Y-axis,  rad.  / 
sec. 

Angular  velocity  about  the  body  Z-axls,  rad.  / 
sec. 

Fan,  induced  lift,  lbs. 

Fan,  induced  moment,  ft.  -lbs. 

Ram  drag  along  the  body  X-axis,  lbs. 

Ram  drag  along  the  body  Y-axis,  lbs. 

Pitch  control  bucket  angle,  deg 

Louver  position,  aft  set,  main  fans,  degrees. 

Louver  position,  forward  set,  main  fans,  degrees. 

Stagger  angle,  main  fan  louvers,  degrees. 

Vane  angle,  main  fan  louvers,  degrees. 

Collective  stagger  stick  position  (positive  for 
increasing  altitude),  ra4i&ns. 

Roll  control  stick  deflection  (positive  for  right 
roll),  radians. 

Collective  vane  stick  position  (positive  for  forward 
acceleration),  radians. 

Yaw  control  stick  deflection  (positive  for  right 
yaw),  radians. 

Pitch  co.itrol  stick  deflection  (positive  for  pitch 
up),  radians. 

Commanded  stagger  angle  per  stick  deflection, 
rad.  /rad. 

Commanded  stagger  angle  per  w,  rad. /ft. /sec. 


K  ^  Commanded  stagger  angle  per  roll  angle,  rad.  / 

rad. 

Kp  Commanded  stagger  angle  per  roll  control 

deflection,  rad.  /rad. 


K  Commanded  vane  angle  per  throttle  position, 

rad.  /rad. 


K 


* 


Commanded  vane  angle  per  yaw  angle,  rad.  /rad. 


Commanded  vane  angle  per  rate  of  yaw,  rad.  / 
rad.  /sec. 


Kg  Commanded  vane  angle  per  yaw  control  deflec¬ 

tion,  rad.  /rad. 

K  v  Commanded  bucket  angle  per  pitch  control  de¬ 

flection,  rad. /rad. 

Kg  Commanded  bucket  angle  per  pitch  angle,  rad./ 

rad. 


Kg  Commanded  bucket  angle  per  rate  of  pitch,  rad./ 

rad.  /sec. 


2. 3  Analysis 

The  forces  and  moments  for  the  hover  flight  regime  were  calcu¬ 
lated  in  the  body  axit,  system,  and  the  body  coordinates  were  referenced 
to  earth  coordinates  by  means  of  a  Yaw,  Roll,  Pitch  set  of  Euler  angles 
and  direction  cosines. 

The  body  axis  equations  of  motion  are  given  below : 

T  T 
xje+  Xr 

u  =  rv  -  qw  + - -jn -  +  X  u  +  J  g  (1) 

U  Z  0 


(Continued) 


w  =  qu-pv  -  (- 


(2) 


T  T 
z.  z 

- — — )  -  T  +  Z  w  +  Zq  +  Zu  +  n  g 

m  nf/m  w  q  u  z  o 


.  I 


I  -  I 


xz  ,22.x  'z  Xmf  Tnf^ nf 


q=—  <r-P>+I -  Pr  +  (TZ  +  T  )  j“"“  +” 

y  y  try 


H  p  Hr 

+ _ - _  +  X  +Mq  +  M  w  +  Mu 

I  I  q  w  u 

y  y 


v  =  pw-ru  +Yv  +  Yr+m  g 
v  r  z  o 


T  T 

I  I  -I  z0  z  „  n 

.  xz  •  v  z  &-  r  H  q 

P“V  (r  +  pq)+-t - qr  +  (-j - ) ”  —t—  +LpP 

x  X  X 


T  T 

I  II  x.  x  h  q 

r  =  y2-  (p-qr)  +-^-  pq  +  (  ~ +  T~  +Nvv  +  Nrr 
z  z  z  z 


(3) 

(4) 

(5) 

(6) 


The  Euler  angular  relationships  used  are  described  in  Section  9. 0 
of  this  report. 

2. 3. 1  Thrust  Functions 

The  fan  thrust  data  was  based  on  the  original  General  Electric 

powerplant  specification.  These  data  were  later  corrected  to  conform 

with  subsequent  data  furnished  by  General  Electric.  For  the  condition 

of  2500  foot  altitude  and  a  standard  hot  day  (p=  .00205), 

T  =  5000-300  (100-%  N J  per  fan. 
o  g  (7) 
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For  small  changes  in  altitude  (<  500  ft. )  the  thrust  was  modified 
by  the  density  ratio  so  that 

Tq  =  [5000-300  (100-%  Ng)  J  [l-.  000026  Ah  j  per  fan  (8) 

The  axial  and  thrust  fan  forces  were  modified  by  vector  and 
stagger  angle  as  below: 

2 

Tz  »  Tq  (1.000189^  )  (1-.  004660  )  (9) 

Tx  =  To  <•  0174-. 0001133  0S)  0V  (10) 

The  nose  fan  thrust  was  assumed  to  be  a  linear  function  of  nose 

fan  bucket  (thrust  reverser)  position,  measured  from  -200  lb.  of  thrust. 

The  nose  fan  thrust  was  assumed  to  vary  linearly  from  -200  to  +  1500  lbs. 

thrust  at  100%  N  for  <5  variation  from  0  to  65° . 
g  P 

The  linear  and  angular  accelerations  due  to  the  thrust  forces  are 
derived  from  the  geometry  of  Figure  2,  and  are  thus  represented  in  the 
body  axis  equation  of  motion. 

2.3.2  Stability  Derivatives 

Since  the  study  was  concerned  primarily  with  small  perturbation 
velocities  around  zero  (  -  50  ft/sec),  the  aerodynamic  characteristic. 3 
of  the  airframe  were  ignored,  and  all  stability  derivatives  were  derived 
from  momentum  considerations  of  the  fan  flow. 

Available  data  was  plotted  to  give  the  linear  variations  of  ram 
drag  with  forward  velocity  as  a  function  of  total  thrust  vector.  Assuming 
the  fan  reaction  forces  to  be  due  to  simple  momentum  change,  mAV,  the 
variation  of  ram  drag  with  velocity  is  simply  the  mass  flow  rate  through 
the  fan.  Based  upon  these  assumptions  the  mass  flow  rate  of  each  main 
fan  is: 


.  214v/Tq  per  fan. 


(11) 


I 


For  the  nose  fan,  using  the  previously  established  relationship 
between  nose  fan  and  main  fan  total  thrust  vector, 

m  r  =  .0545  J  T 

nf  v  o  (12) 

The  mass  flow  rate  was  assumed  constant  at  a  nominal  thrust  cf 
4850  lbs/fan  which  was  approximately  that  required  to  hover,  so  that: 

m  ,  =  3. 86  slugs/sec.  (13) 

ni 


m  ,  =  15. 15  slugs/sec. 
mf 


The  forces  due  to  momentum  changes  can  be  expressed  as: 


Fj  =  -  (Vj  +  w  X  (15) 

Where  SXRt  is  the  velocity  due  to  rotation  u,  and  R^  ■  x^t  +  +  z^k. 

In  body  axes,  then 

—  •  r  A  A  A  _  _  , 


Fj  =  -  Jul  +  vj  +  wk  +  wXBj  J 


t  t  £ 


—  _  A  /s 

Now,  wXRj  =  p  q  r  =i  (qz^ry^  +  j  (rx^-pz^) 


xi  yi  zi 


+  k  (py1-qxJ) 


in  our  case,  z^  =  0,  so  eq.  (16)  becomes,  when  expanded, 
Fj  =  mi  [  f  (u-ryp  +  f  (v  +  rx^  +ic  (w  +  py^qx^  j 
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For  our  system,  the  fan  coordinates  are: 


= 

15.45' 

xmfr 

= 

-1.0 

= 

-1.0 

ynf 

= 

0 

Vfr 

= 

+5 

yml  l 

= 

-5 

So  that,  if  we  define  the  dimensional  derivative  as  acceleration/velocity; 


m  2m 

X  =  Y  =  Z  =-  - — *)  =  118/sec. 

u  v  w  m 


(19) 


m  fx  Z  . 
nf  nf  +  m  -x  , 

Z  =  -Y  = _ I _ 101/sec. 

q  r 


(20) 


m 


The  moments  due  to  these  forces  can  be  expressed  as  R.  X  F.. 
Thus,  taking  the  cross  product  of  R^  and  equa.  (18),  F^,  we  get 


t 


1 


A 

k 


R.  XF.=m 
i  i  i 


*i 


\ 


(21) 


(u-ryi)  (v+r  x.)  (w  Hpy.-qx^ J 


and  expanding, 

¥V  -t 


m 


1 


|yt  ( W  +  py4  -  qx^  -  z{  (v+TXj) 


+  J  \\  (u-ryj)  -  x.  (w  +  py^qx^ 


A 

+  k 


and  letting  =  0 


*t  (v+rxt)  -  y{  (u-ry^ 


M  =  -m 


f  (wyt  +  py^-qx^)  +  f  (-wxj  +  qxi2-px.yi) 


a  2  2 

+  k  (vxA  +  rxA  -uyi  +  ryt  ) 


(22) 


(23) 


and  we  can  calculate  the  dimensional  derivatives  for  rotation 


L  = 


2mmfXmf 
I 


=  -.0427/sec. 


x 

1nfXnf 


m  jc  .  .  2m 


M  *. 

w 


+  x  . 

mf  mf 


=  00172/ft.  sec. 


(24) 


(25) 


^/rf  +  2*mfx2mf 


M  = 

q 


=  -.0562/sec. 


(26) 


mnfx„f  ,  2mmfxmf 

N  =  -  nfnf+  mf  mf  =  -.00155/ft. /sec. 
v  ; 


(27) 


•  2  « •  ,  2  A  2  . 
m  jc  -  +  2m  .  (x  _  +  y  ,) 

VT  nf  nf  mf  mf  mf  --A_  /oov 

N  =  -  ■■  =  -.  0906/ft.  sec.  (28) 

r 

I 
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and  because  of  symmetry  in  the  x  -  y  plane: 


X  =  Z  =  L  =  L  =  M  =  N  =0. 
r  p  w  q  p  u 

The  only  additional  stability  derivatives  considered  were  those  due 
to  induced  lift  and  those  due  to  control  motion.  The  induced  lift  and  mo¬ 
ment  derivatives  were  taken  as: 


Z  =  -.275  =  _>0954/sec>  at  4850  lbs>  thr<  (29) 

U  M 

1.078  /2T 

M  =  - - - -  =  .00637/ft.  sec.  at  4850  lbs.  thr.  (30) 

u  I 

y 

The  control  derivatives  were  continuously  computed  as  shown  in 
equations  (9)  and  (10).  For  small  perturbations  around  hover,  and 
the  thrust  per  fan  of  4850  lbs. 


57.3  X  .00932  T 


o 


/? 


m 


r.  9. 1  l't/sec" 


(31) 


I. 


/3 


:  diff 


37.3  X.00MSV  2 

- — -  =  3.71/sec 

I 

x 


(32) 


57.3  (.0174-.  0001133  /?  )  m  ,2 

M  .  s'  2v  .  T  =  2.4/sec 

ii  - *  mf  o 

v  I 


(33) 


This  is  calculated  at  12°  P  on  each  fan,  which  is  required  to 


hover  at  4850  lbs.  thrust. 


Similarly  - 


(.  0174-.  000113  0  )  2T  o 

X„  =  57.3 .  8 - 0  =  31.3  ft/sec  (34) 


m 


„  57.3  (26.2)  xo  ,  2 

Z6B  =  - - -  ^ -  =  -4.17  ft/sec 


m 


100% 


(35) 


M6 


B 


57.3(26.2)x  T  « 

-S^-2-=  1.11/sec 


I  T 

y  o 


100% 


(36) 


2.3.3  Control  Command  Equations 

The  command  equations  comprise  the  total  control  surface  pos¬ 
ition  commanded  as  a  function  of  the  inputs  from  the  pilot  and  from  the 
SA  system  sensors. 

The  SA  system  utilizes  a  three  axis  rate  gyro  package  and  asso¬ 
ciated  networks  to  provide  "pseudo  position"  and  rate  feedback  signals. 


The  SA  compensation  networks  have  the  transfer  function, 

(1  +  RrS)  (37) 


60  _  =  K 

61 


1  +TS 


And  thus  the  total  control  command  equations  can  be  written: 


(l+R/S) 

6  =  K  v+  K  _ i _ q 

Bc  9  1  +  T  S 


(38) 


0 


-  K„<r-  (1  +  Vs*  p  +  K„r 


1  +  TS 


(39) 


(40) 


XV  ^  W 

*v  =  KA^f  1  +  T  {J  )  r  "  Kj;^ 

C 

r 

V  ■  KrT+  K*  p  -  Krr  (41) 

Ci 

*v  ~  KAX~  K»  <IV+4^s‘S>  r  »  Ktt  <42> 


Control  motions  commanded  by  the  SA  system  were  limited  to 
25%  of  the  total  authority,  around  trim.  These  limits  were: 

a<5_  =  ±8.125  deg. 

BLim 

=  -  3. 5  deg. 

SLim 

=  -  3  deg. 

Lim 

Figures  3,  4,  and  5  are  the  computer  diagrams  used  to  mech¬ 
anize  the  hovering  equation  of  motion. 


22 


■*'  P  >  l>h 


I 


'J' 

L  ■  Ur  | 


Figure  3  Computer  Diagram 
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3.0  LONGITUDINAL  THREE-DEGREE-OF- FREEDOM 
TRANSITION  SIMULATION 


3.1  Discussion 


The  basis  for  the  longitudinal  transition  work  described  in  this 
section  was  the  2500  foot,  hot  day,  flight  condition;  an  aircraft  with  the 
eg  at  station  246;  and  a  weight  of  9200  pounds.  The  analog  computer 
simulation  work  Involved  the  use  of  a  fixed  base  oscilloscope  display 
cockpit  for  providing  pilot  inputs  to  the  simulation,  as  well  as  for  obtain¬ 
ing  pilot  evaluation  of  system  operation. 


This  simulation  included  investigation  or  development  of  the 
following  items: 

1.  Altitude  Control  and  Transition  Procedures.  Pitch  control 
with  and  without  stabilization  system 

2.  Single  engine  recovery  and  speed-altitude  envelope  for  10  fps 
sink  rate 


3.  Selection  of  vector  rate  for  pilot  command  during  transition 

4.  Trimmed  flight  at  various  a,  (3^  and 

5.  Methods  of  controlling  horizontal  tail  trim  during  transition 

6.  Nose  fan  phase-out  with  /?  ,  evaluation,  and  nose-fan  and  tail 

*  y 

malfunctions  during  transition 

7 .  Variation  of  d  Cp  /ddp  on  handling  qualities 

8.  Pitch  auto  stabilization  authority,  gains,  and  concept  in  the 
transition  regime 

9.  Hand  calculation  dynamic  check  of  simulation 
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The  mechanics  of  the  simulation!  as  well  as  the  limitation  of  the 
simulation  are  included  in  the  following  Section  3. 15.  (Shown  are  approx¬ 
imation.  data,  data  fits  (separation  of  data),  limitations,  etc.) 

3.2  Accelerating  Transition 

The  accelerating  transitions  require  that  the  control  stick  be  dis¬ 
placed  to  a  nose-down  position  during  the  entire  sequence,  which  in  turn 
means  the  maneuvering  gain  switch  on  the  stick  will  be  energized.  The 
maneuvering  gain  in  the  pitch  mode  must  thus  be  chosen  to  afford  damp¬ 
ing  according  to  the  pilot's  wishes.  In  the  transition  simulation,  ener¬ 
gizing  the  maneuvering  switch  cuts  the  pitch  channel  gain  to  zero.  The 
transition  offered  no  control  problems  in  the  absence  of  auto-stab  signals . 
so  the  results  should  be  conservative.  Figure  6  is  a  plot  of  %  stick  po¬ 
sition  versus  velocity  for  trimmed  velocities  at  various  0^* 

A  fast  accelerating  transition  required  more  forward  stick  than 
a  slow  accelerating  transition,  because  in  the  fast  transition  higher  vector 
angles  are  obtained  at  lower  speeds,  which  means  the  tail  contribution 
opposing  the  pitch-up  tendency  is  less.  The  "fast"  transition  required 
twice  the  forward  stick  for  a  trimmed  transition. 

There  are  two  limits  placed  upon  the  maximum  allowable  vector 
rate  during  the  accelerating  transition.  They  are:  (1)  the  fact  that  the 
higher  the  vector  rate,  the  more  forward  stick  is  required  to  hold  the 
nose  down,  and  (2)  the  fact  that  a  higher  vector  rate  results  in  a  higher 
thrust  level  being  required  to  hold  lift  equals  weight.  If  the  gas  gener¬ 
ators  are  already  at  100%  and  the  louvers  are  already  at  maximum 
lift,  the  vector  rate  must  be  slow  enough  so  the  induced  and  aerodynamic 
lift  components  increase  sufficiently  to  compensate  for  the  loss  of  lift  due 
to  thrust  vectoring.  In  general,  the  maximum  vector  rate  limit  due  to 
insufficient  lift  occurs  before  the  limit  due  to  insufficient  longitudinal  con¬ 
trol  power. 

To  determine  the  maximum  allowable  vector  rate  during  an  accel¬ 
erating  transition,  a  high  gain  autopilot  was  synthesized  and  connected  in¬ 
to  the  pitch  mode.  Pitch  angle  was  fed  to  the  control  stick  input,  and  rate- 
of-climb  was  fed  to  collective  stagger.  Ngwas  set  at  100%.  Collective 
vector  was  set  to  operate  at  a  preset  rate.  The  computer  was  then  turned 
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on,  and  several  transitions  were  made  to  determine  the  maximum  vector 
rate  which  required  no  less  than  13°  collective  stagger  and  resulted  in 
no  change  of  altitude  or  pitch  attitude.  As  simulated,  this  maximum  vec¬ 
tor  rate  was  found  to  be  2°  per  second. 

Other  considerations  dictated  a  faster  rate  than  this,  and  3°  per 
second  was  chosen.  This  means  that  the  pilot  must  "beep"  the  vector  but¬ 
ton  to  attain  an  average  vector  rate  of  around  2°  per  second,  maximum. 

3. 2. 1  Preferred  Technique  for  Accelerating  Transition 

The  preferred  accelerating  transition  technique  was  developed 
to  be  as  follows: 

1.  Set  stagger  at  maximum  lift. 

2.  Increase  throttle  to  100%  N  . 

g 

3.  After  lift-off  immediately  begin  vectoring  until  rate  of  climb 
goes  to  zero  on  IVSI. 

4.  Control  rate-of-climb  with  reference  to  the  IVSI  ( Instantaneous 
Vertical  Speed  Indicator )  by  beeping  vector  until  rate-of-climb 
becomes  negative.  As  rate-of-climb  goes  to  zero,  again  resume 
beeping. 

Using  this  procedure,  accelerating  transitions  can  be  consistently  made 
within  a  -  10  ft.  altitude  band. 

3. 3  Decelerating  Transition 

The  decelerating  transition  requires  a  different  approach.  Dur¬ 
ing  this  maneuver,  the  aerodynamic  lift  is  appreciable,  and  it  is  indeed 
necessary  to  reduce  thrust  to  slow  down.  Reduction  of  gas  generator 
RPM  also  reduces  the  nose-up  pitching  moment  caused  by  the  fans  at  in¬ 
termediate  vector  angles,  and  thus  allows  more  rapid  removal  of  collect¬ 
ive  vector.  The  decelerating  transition  can  be  accomplished  at  the  max¬ 
imum  3° /second  vector  rate. 

The  preferred  decelerating  transition  procedure  was  developed 
to  be  as  follows: 

1.  Reduce  speed  to  conversion  speed  at  90%  N^. 
la.  Set  lift  control  at  maximum  lift. 
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2.  Convert. 

3.  Vector  towards  j3  =0,  until  a  rate-of-climb  shows  on  the  IVSI, 

v 

4.  Cease  vectoring  until  a  sink  rate  appears  on  the  IVSI,  and  then 
resume  vectoring.  Repeat  this  operation  as  often  as  required. 

5.  At  40-50  knots  IAS  begin  applying  power  to  keep  sink  rate  zero. 
The  vector  angle  will  be  around  zero  at  this  time. 

6.  Hover  or  land. 

The  same  1 10  ft.  altitude  corridor  is  possible  with  this  transition  also. 

These  techniques  result  in  the  fastest  transition  possible  with 
zero  altitude  loss  or  gain.  Transition  may  be  made  at  as  low  a  rate  as 
desired,  and  lower  rate  transitions  require  less  forward  stick. 

The  fastest  accelerating  transition  requires  80%  forward  stick  at 
intermediate  vector  angles,  whereas  a  slow  transition  which  does  not  de¬ 
viate  appreciably  from  thrust=drag  can  be  accomplished  with  40%  forward 
stick. 


The  decelerating  transition  requires  less  forward  stick  because  it 
is  accomplished  mostly  at  higher  speeds  where  the  tail  and  elevator  effec¬ 
tiveness  are  appreciable. 

True  hovering  (-  360°O!)  was  not  possible  in  this  simulation  due  to 
the  method  used  of  resolving  body  rates  into  angle  of  attack.  Due  to  the 
extensive  hovering  simulation  already  conducted,  it  was  felt  that  rearward 
flight  was  not  important  to  the  transition  investigations.  The  system  used 
allowed  true  vertical  flight,  with  a  being  limited  to  ±  1400  . 

3.4  One- Engine  Out  Operation 

For  evaluation  of  single-engine  fan  powered  flight,  the  fan  thrust 
data  was  modified  to  give  59%  of  the  two-engine  main  fan  thrust  and  56% 
of  the  two-engine  nose  fan  thrust. 

A  speed- altitude  envelope  was  determined  by  flying  at  a  given 
speed  and  zero  sink  rate,  and  then  reducing  fan  thrust  by  means  of  a 
switch  to  the  single  engine  value,  taking  into  account  the  fan  dynamics. 


An  indicated  airspeed  of  45  knots  or  above  allows  10  fps  sink  rate  or  less 
to  be  established  and  held.  The  criterion  of  a  successful  recovery  was 
to  determine  how  much  altitude  must  be  lost  to  reduce  the  sink  rate  to  ten 
feet  per  second.  The  technique  of  recc  *ery  varies  with  the  initial  altitude. 
For  the  high-altitude  recovery,  assuming  low  vector  angle  and  low  veloc¬ 
ity,  the  recovery  technique  involves  immediate  increase  of  the  remaining 
Ng  to  100%,  dropping  of  the  nose,  and  immediate  initiation  of  vectoring 
toward  45°  vector  angle.  From  ';he  high  altitude  one-engine  out  recovery 
aspect,  it  is  evident  that  the  highest  vector  rate  possible  is  desirable. 

The  low  altitude  recovery  in  the  hovering  regime  involves  merely 

applying  maximum  power  to  the  remaining  engine  and  holding  on.  At  30 

knots  or  so  IAS,  it  is  desirable  to  vector  towards  45°  0  and  attempt  to 
,  v 

gain  velocity. 

A  plot  of  the  one-engine  out  speed  altitude  envelope  for  10  fps 
sink  rate  is  shown  in  Figure  7. 

3. 5  Choice  of  Pilot  -  Commanded  Collective  Vector  Rate  During  Fan- 
Powered  Flight 

As  has  been  discussed,  two  conflicting  requirements  have  affected 
the  choice  of  collective  vector  rate  during  the  transition.  They  are: 

1.  The  rate  of  vectoring  during  the  accelerating  transition  is 
limited  by  the  requirement  to  be  able  to  perform  zero  alti¬ 
tude  change  transitions.  This  limit  is  2°  per  second. 

2.  A  higher  rate  of  vectoring  is  useable  during  the  decelerating 
transition,  and  an  even  higher  rate  of  vectoring  is  desirable 
to  accomplish  a  high-altitude  one-engine  recovery. 

It  was  decided  that  the  pilot  should  have  a  capacity  of  vectoring  in  excess 
of  the  maximum  allowable  for  the  accelerating  transition,  but  that  the  ex¬ 
cess  should  be  limited  so  that  all  inadvertent  steady  vectoring  to  maximum 
vector  would  not  be  catastrophic. 

3.6  Computer  Derived  Trim  Conditions 

During  the  longitudinal  transition  simulation,  a  number  of  trim 
conditions  were  arrived  at  by  manually  "flying"  the  simulation  to  a  given 
indicated  air-speed  at  a  given  vector  angle,  and  then  going  to  "hold"  on 
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the  computer  and  reading  the  various  flight  parameters  out.  It  should  be 
noted  that  for  this  simulation  sin  6  was  defined  equal  to  0  in  radians,  so 
high  angles  of  pitch  (and  a  )  resulted  in  errors  in  the  trim  condition.  The 
data  is  reasonable  up  to  ±  io°  angle  of  attack.  The  results  of  this  work 
are  plotted  in  Figure  8. 

3.7  Tail  Trim  During  Transition 

Several  schemes  to  position  the  horizontal  tail  during  the  transi¬ 
tion  were  investigated. 

For  the  accelerating  transition,  a  considerable  nose-up  moment 
from  the  fans  must  be  balanced  by  control  inputs.  Figure  9  is  a  com¬ 
puter-derived  plot  of  required  tail  incidence  versus  velocity  for  trimmed 
fan-powered  flight. 

It  would  appear  that  a  simple  program  of  tail  incidence  versus  0  y 
would  suffice  for  transition  operation  (untrimmed  fan-powered  flight). 

This  is  not  true  however,  since  the  vector  angle-velocity  relationship  is 
different  for  accelerating  and  decelerating  transitions.  Furthermore,  the 
decelerating  transition  occurs  at  a  different  gas  generator  RPM  than  the 
accelerating  transition.  The  tail  incidence  turns  out  to  be  a  function  of 
velocity  and  power  setting,  and  since  transitions  are  made  under  varying 
flight  conditions,  an  open- loop  program  of  tail  versus  /?v  is  unsatisfac¬ 
tory.  Many  open-loop  tail  programs  were  tested,  and  they  universally 
were  wrong  for  one  or  the  other  of  the  transitions.  A  program  of  tail 
incidence  which  was  a  function  of  vector  angle  and  power  setting  worked 
reasonably  well  at  the  2500  foot  hot  day  condition,  but  would  be  wrong  for 
S.  L.  std.  day  conditions. 

For  the  above  reasons,  programmed  open-loop  tail  incidence  was 
discarded.  Some  form  of  closed  loop  (pilot  input)  scheme  was  indicated, 
and  the  following  were  tested: 

1.  Pilot  manual  Beeper  control. 

2.  Integral  control  based  on  returning  the  elevator  position  to 
zero  in  the  steady  state. 

3.  Automatic  beeper  control  using  elevator  position  by  means 
of  control  stick  position-actuated  switches.  Using  this  sys¬ 
tem,  the  tail  moves  at  a  fixed  rate  if  the  control  stick  is  dis¬ 
placed  more  than  a  predetermined  amount. 
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4.  Automatic  beeper  control  using  a  combination  of  elevator 
position  sensing  and  stability  augmentation  system  output 
sensing.  This  system  commands  a  fixed  tail  rate  if  the 
sum  of  control  stick  plus  auto  stab  output  ($  )  exceeds 

a  certain  predetermined  value.  a 

The  Cooper  pilot  rating  was  used  to  classify  these  various  systems.  Dur¬ 
ing  these  runs,  the  pilot  made  no  specific  maneuvers,  but  just  gave  opin¬ 
ions  on  his  feelings  toward  the  handling  qualities.  The  Ryan  Racer,  a 
fixed  base  cockpit,  was  used  for  these  tests. 

In  general,  the  results  of  these  tests  were  as  follows: 

1.  The  manual  beeper  system  had  a  best  rating  of  4  obtained  with 
a  3°  /sec  tail  rate.  Other  tail  rates  had  poorer  ratings. 

2.  Integral  control  had  a  best  rating  of  3. 

3.  Elevator  alone  beeping  the  tail  had  a  best  rating  of  3-4 

with  a  1°  6  dead  band  and  2°  /sec  tail  rate, 
e 

4.  Elevator  alone  with  a  two  rate  system  and  two  dead  bands, 

a  small  one  for  the  small  rate  and  a  large  one  for  the  larger 
rate,  had  ratings  of  5  or  worse. 

5.  A  combination  of  6e  and  <5p  commanding  a  given  tail  rate 
when  the  sum  of  the  two  exceeded  a  given  dead  band,  had  a 
best  rate  of  3. 

As  a  result  of  these  tests,  the  method  used  in  test  4  was  discarded.  Since 
the  method  of  test  2  would  mean  a  different  tail  trim  system,  and  the  Cooper 
rating  was  not  appreciably  higher  than  tests  1,  3,  and  5,  the  integral  con¬ 
trol  was  discarded. 

The  method  of  test  5  gave  good  results,  but  required  the  mixing  of 
electrical  (auto  stab  output)  and  mechanical  (control  stick  position)  infor¬ 
mation.  It  was  not  discarded,  but  just  put  in  abeyance. 

Hie  manual  beeper  system  and  the  control  stick  switch  beeper  sys¬ 
tem  had  about  the  same  rating,  but  method  3  had  certain  advantages.  They 
are: 
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1.  Due  to  the  requirement  of  forward  stick  displacement  during 
transition,  the  horizontal  tail  is  automatically  positioned  at 
the  full  up  position  by  the  time  it  begins  to  be  effective. 

'  2.  The  requirement  of  pilot  manual  beeping  would  detract  from 
the  ability  of  the  pilot  to  hold  zero  rate  of  climb  by  selective 
vector  beeping.  He  only  has  one  thumb  on  the  right  hand. 

The  control  stick  beeping  allows  the  pilot  to  devote  nearly  full 
time  to  his  transition  requirements. 

Using  either  system  1  or  system  2,  it  is  possible  to  trim  the  airplane  with 
stick  centered  so  that  full  auto  stab  nose  fan  authority  is  being  commanded. 
If  the  stab  system  is  thus  saturated,  the  auto  stab  system  is  effectively  im¬ 
perative  until  the  hardcover  bleeds  away.  Due  to  the  high  system  gain  pos¬ 
sible,  this  effect  could  easily  occur.  On  the  other  hand,  this  problem  only 
presents  itself  at  the  higher  transition  speeds,  whr  re  auto  stab  operation 
is  not  really  required. 

System  5  eliminates  this  effect,  but  until  evidence  to  the  contrary 
is  available,  this  effect  will  be  considered  unimportant. 

System  3  has  been  deemed  the  most  attractive  on  the  combined  basis 
of  simplicity,  reliability,  and  operational  effectiveness. 

3.8  Nose  Fan  Phase-Out  With  P  v 


As  vector  angle  is  increased,  the  pilot  stick  input  to  the  nose-fan 
thrust  reverser  doors  is  phased  out,  and  the  nose  fan  thrust  level  goes  to 
zero,  although  the  auto  stab  inputs  still  have  the  original  25%  of  stick  au¬ 
thority.  Transitions  were  accomplished  with  various  programs  of  pitch 
fan  phase-out  with  vector  angle,  and  no  difference  could  be  noted  in  air¬ 
craft  response.  For  this  reason,  the  pitch  fan  phase-out  program  incor¬ 
porated  into  the  airplane  was  deemed  satisfactory. 

3.  9  Tail  Incidence  Actuator  Malfunction  During  Fan- Powered  Flight 

If  the  tail  actuator  fails  hard  over  in  the  up-trim  position,  nose 
down  control  power  runs  out  at  20°  (3  v  (55  Knots  EAS). 

If  the  actuator  fails  in  the  nose-down  trim  position,  the  transition 
can  be  successfully  completed,  but  considerable  back  stick  is  required. 


The  minimum  tail  incidence  allowable  for  a  successful  transition 
is  from  5°  to  7°  up  . 

If  the  tail  fails  full  up  at  higher  than  conversion  speed,  sufficient 
back  stick  might  not  he  available. 

3. 10  Nose  Fan  Door  Failure 


In  the  case  of  a  nose  fan  door  failure,  the  vehicle,  of  course,  can¬ 
not  be  hovered.  For  a  hard-over  in  the  maximum  thrust  direction,  stick 
control  runs  out  at  75  Knots  EAS. 

If  the  doors  fail  in  a  closed  position  while  the  fan  is  running,  we 
are  indeed  in  deep  trouble. 

A  failure  resulting  in  zero  nose  fan  thrust  results  in  all  pitch  con¬ 
trol  resulting  from  the  elevators  and  tail;  30  Knots  would  be  considered 
a  minimum  speed  for  this  condition,  and  control  is  very  mushy. 

3.11  Variation  of  $C  on  Handling  Qualities 

<p 

This  coefficient  causes  the  nose  fan  RPM  to  vary  with  nose  fan 
thrust  reverser  door  position.  The  value  of  this  coefficient  was  increased 
by  a  factor  of  10  with  no  apparent  effect. 

3. 12  Pitch  Channel  Auto  Stab  Authority.  Gains  and  Concept  in  the 

Transition  Regime 

Of  all  the  work  to  date,  the  transition  regime  seemed  to  be  the 
least  critical,  at  least  in  the  longitudinal  mode  as  far  as  the  auto  stab 
system  is  concerned.  Most  of  the  transitions  were  accomplished  with 
maneuvering  gain.  Since  the  stick  is  displaced  outside  the  holding  gain 
limits  for  practically  the  entire  transition,  maneuvering  gain  was  set 
to  zero  for  all  transition  work.  At  transition  flight  trim  speeds,  the 
dynamic  response  of  the  airplane  alone  and  airplane  and  auto  stab  was 
recorded.  The  autopilot  reduced  the  characteristic  divergent  character 
of  the  airplane  response  to  a  considerable  degree. 

3. 13  Dynamic  Check  of  Simulation 

As  with  all  analog  simulation  work,  some  sort  of  independent 
dynamic  response  check  is  required  if  the  maximum  confidence  of  valid 
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results  is  to  be  obtained. 


Due  to  the  difficulty  of  calculating  this  check,  a  single  case. 

35°  v.  100  %  RPM  at  angle  of  attack  equals  zero  was  picked  for 
calculation.  The  calculations  were  made  by  a  small  perturbation  anal¬ 
ysis  around  trim  points.  For  this  condition  e  / 6  p  was  derived  to  be: 


-1.32  (S+.227)  (S+,628) 

(S-.219)  (S+.45)  2  o 

(S*  +  2(.  614)  (1.32)  S  +  1.32  ) 


The  transient  response  of  this  system  to  an  impulse  of  6  p  was 
derived  to  be  as  follows;  for  a  unit  6  p  impulse 


r  .219t  -.  45t 

6  (t)  =  -1.32  l.2639e  +  .0492e 

+  .752e"*81t  Sin  (1.04t  -  24.6°)J 

A  plot  of  this  response  is  shown  in  Figure  10  and  a  plot  of  the  computer 
derived  response  for  the  same  condition  is  shown  in  Figure  11. 


STICK  AFT 
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Figure  10  Dynamic  Check  Hand- Calculation  Response  for  impulse  o i  6  p,  ior  Small  Pert. 
Transfer  Function 


TIME 

Figure  11  Plot  of  Analog  Computer  Record  for  Dynamic  Check  Initial  Conditions 
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3.14 


List  of  Symbols 


i 


\ 


u 

w 

e 

A 


A 


F 

x 


F 

z 


M 


m 


g 

I 

y 

L 


Velocity  along  vehiole  X  axis,  ft/sec. 
Velooity  along  vehiole  Z  axis,  ft/sec. 
Pitch  angle,  radians 

Time  derivative  of  arbitrary  quantity,  A. 
Second  time  derivative  of  A 
Foroe  along  vehicle  X  axis,  lbs. 

Foroe  along  vehiole  Z  axis,  lbs. 

Pitching  moment,  ft.  lbs. 

Mass,  slugs. 

32.2  ft/sec2 

2 

Pitch  moment  of  inertia,  slug  ft 

Lift,  perpendicular  to  velocity  vector/lbs. 

% 

Total  velocity,  ft/sec. 


Slipstream  dynamic  pressure,  q  +  Tooo/, 


P/2VT  •  PSP 


Density  of  air,  slugs/cu.ft. 

2 

Total  wing  fan  area,  ft 


Slipstream  normal  force  coefficient  - 

q8a 


PSF. 
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N 


N 


Normal  force,  -F  ,  lb. 

z 


Tail  moment  arm,  feet. 


s 


CVT  at  a  =  0 
N 


PF 


Pitch  fan  area,  ft2;  and  Axial  force,  lb. 


~f 

a 


0 


Main  fan  diameter,  ft 


Angle  of  attack,  degrees  or  radians,  as  specified, 
Collective  vector  angle,  degrees. 


0 


8 


R 


Collective  stagger  angle,  degrees. 


1-T 


s 


ooo 


000/ 


AF 


q  +  T 


000/af 


Wing  fan  thrust  at  0^=  0g=VT=O 


NF 


m 


M 


PF 


Axial  Force  =  -Fx,  lbs. 


Slipstream  axial  force  coefficient 


q  Af 

Nose  fan  ram  drag,  considered  as  a  body  axis 
term  here,  lbs. 

M 

Slipstream  pitching  moment  coefficient  "g 


q  A{Df 

Pitch  fan  pitching  moment,  MQ  +  Npj^pp  ft 

PF 


lbs. 


a.'5 


Pitch  fan  pitching  moment  about  fan  hub  ft.  lbs. 


e 


«e 


**6 


Ma 


N 

g 

6f 

N 

T 

g 


Downwash  angle  at  horizontal  tail,  degrees. 

Tail  incidence  angle  w/respect  to  body  X  axis, 
degrees. 

Elevator  deflection,  degrees. 

2 

Horizontal  tail  area,  ft 


Tail  moment  arm,  ft. 
Damping  moment  due  to 
Damping  moment  due  to 
Fan  power  coefficient—. 
Fan  power. 

Gas  generator  RPM,  % 


6  ,  ft.  lbs/rad.  /sec. 
a,  ft. lbs. /rad. /sec. 


Flap  deflection,  degrees 
%N 

i— !_ 

100 


Gross  thrust 


3. 15  Analysis 


The  equations  of  motion  used  for  the  longitudinal  three-degree- 
of-freedom  transition  simulation  are  shown  below.  The  equations  are 
in  body  axis  form,  and  are  related  to  the  inertial  axes  by  the  pitch 
angle,  6  . 

Z  F 


ti  =  x  -  gsinfl  -  qw  (1) 


m 

m 

ZF 

w  =  _JL 

+  geos  6  +  qu 

(2) 

m 

m 

..  ZM 

9  =  1 

y 

(3) 

Where: 

2F  =  -L 
x  tail 

a 

sin  (  2  "  €  0  ) 

-  apf  -  ca“af 

(4) 

Z  F  =  -L 
z  tail 

cos  (*f  _  c0  ) 

-  Cx.SqSA  -  N__ 

N  F  PF 

(5) 

Z  M  -  MT1T, 
PF 

s  s 

+  C  q  AD  - 
m  f  f 

99qu  -  87  5q 

(6) 

•Vt 


where  -99qu  is  equal  to  the  tail  pitching  moment  due  to  pitch  rate,  and 
-87 5q  is  the  damping  moment  due  to  the  fans.  The  drag  of  the  tail  is  ig¬ 
nored  as  a  pitching  moment  contribution. 


3.15.1  NORMAL  FORCE 


FXT  =  -  F„ 
N  Z 


C  -  C  +AC 
N  N  N 

o  ct 


+  ACN  8  +acn/ 

0  p  s 

^  V 


(7) 
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Where 


CN  8  =  CN  *  Q  less  the  ACNB  due  to  45°  flap  (Figure  12,  Figure  13) 
o  a~ 


A  C.T  =  Change  in  normal  force  coeff.  due  to  angle  of  attack. 
Na 


V  '  -CN08  -  K0 

0  0  v 

v  L  L 


(Figure  14) 


s  s 

AC  =  -C  K 

N  N_  a  T 

0  O  p  sL 

s 


(Figure  15) 


NpF  -  Normal  force  due  to  the  pitch  fan 


=  3  R 


3.15.2  AXIAL  FORCE 


F  =  -F 
A  x 


CA  =CA  +ACa  a  +ACa  £v  +  ACa  0s 

A  A  Att  A  v  Aa  8 

o  a  B  p8 

yv 


_  8  _  S 

—  C^  at  ot  *  0 
o 


8  8 
AC,  =  change  in  CA  due  to  a 

A  of  A 


(Figure  12,  Figure  13) 
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(Figure  14) 


(13) 


:N  “  > 

o  V  . 


AC 


s 


=  -  C 


P 


N8  V 

o  hb  , 


(Figure  15)  (14) 


DpF  =  Nose  Fan  Ram  Drag  (Figure  16) 

3.15.3  PITCHING  MOMENT  (At  c.g.  Sta.  246) 

M  =  C  SqSA  D,  +  M  +  +  M.q  +  M .  at  (15) 

m  F  F  t  PF  9  a 


C  8  =  C  8  +  AC  8<*+  AC  8  p  +  AC  8  0 

m  m0  “«  mfi  V  8 

v  s 


(16) 


0 


(Figure  17,  18) 


s 

AC  =  Change  in  moment  coeff.  due  to  a  (Figure  17,  18) 
ma 
s 

AC  =  vector  pitching  effectiveness  (Figure  19) 
m 


M  =  Tail  pitching  moment 

■  — — — 

Mpf,  =  Nose  fan  pitching  moment. 

M  •  =  Damping  moment  due  to  fans  and  horizontal  tail. 

U 

=  Moment  due  to  rate  of  change  of  a . 


3.15,4  Force  and  Moment  Equations  for  the  Horizontal  Tall 

For  the  transition  simulation,  the  drag  of  the  tall  was  assumed 

negligible,  and  the  only  contribution  to  pitching  moment  from  the  tail 

was  the  tail  lift.  This  also  held  true  for  the  tail  contributions  to  F 

and  F  .  X 

z 


With  this  in  mind,  the  equation  for  lift  and  pitching  moment 
from  the  tail  are: 


do 

d<  =  .5 

do 

€q  =  f  (T,8)  (Figure  21) 

il  incidence 

=  . 027/deg. 
dfi 

e 

V±25° 

S  *  50. 66  ft2 
i  =  20.6  ft. 

V 
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I 

I 


I 


The  tail  lift  force  is  in  tall  stability  axes,  which  is  removed 
fiom  the  body  axis  system  by  the  angle  a  which  is  equal  to 

de 

Since  =  .  5,  the  angle  of  attack  of  the  tail  is  equal  to 


We  can  thus  write. 

F  =  -L.  cos 

z  tail  t 

(f -*o) 

(19) 

F  =  +  L  sin 

x  tail  t 

(*-.) 

(20) 

M  =  -20.66  L 

(21) 

The  damping  terms  in  pitch,  M  •  and  M  •  are  given  below: 

V  U 

=  -66V.J,  ft  lbs/rad/sec. 

=  -33Vt  ft  lbs/rad/sec. 

Where  is  total  velocity  in  feet/second. 

•  • 

We  assumed  a  »  6  ,  so  we  could  multiply  the  sum  of  M  •  and 
by  V^,q  to  get  the  damping  moment. 

3.15.5  Calculation  of  Thrust  Coefficient  for  the  Fans  and  Forces 
and  Moments  for  the  Nose  Fan 


Main  Fan  Thrust  Coefficient,  T  8. 
_ c 


(Continued) 


The  main  fan  thrust  coefficient, 


000/a, 


_V_  +  T 

P  T  000/ 

2  AF 


T  8,  Is  defined: 
c 


(22) 


Where  T000  is  defined  as  the  thrust  which  would  occur  at  a  given  fan 
RPM  if  0  ,  p  and  Velocity  were  suddenly  all  made  equal  to  zero. 

V 

In  slipstream  notation,  q8  is  defined  as  the  free  stream  dynamic  pres¬ 
sure  plus  the  disc  loading. 

A 

It  can  be  easily  shown  that  ■  ■  1 ..  t  8 

s  c 

q 

8 

1  -  T  has  been  called  R  for  simulation  purposes.  This  is  the  ratio 
c  q 

of  free- stream  to  slipstream  dynamic  pressures. 


The  curves  presented  in  Figures  22,  23,  and  24  are  required  in  the 
computation  of  Tc®. 


g 

The  equation  governing  C  and  T  is: 

p  000 


(23) 


8  I 

C  is  a  function  of  T 
P  0 

P„  is  a  function  of  N  , 
F  g 


,  p  and  p  ,  as  shown  in  Figures  22  and  23. 

V  s 

as  is  shown  in  Figure  24. 


We  can  solve  equation  (23)  for  T000/A  ^ 

2/3 


l000 


PFP 


1/2 


C  8  A 
P  t 


(24) 


We  note  also  in  Figure  24  that  gross  thrust,  T  ,  is  equal  to: 

o 


p  p 

_ E _ 


1/2  \  2/3 


cY 

p  f 


i.  c. 


By  assuming  variations  of  T000/  around  some  initial  value,  we  can 

Af 

write  T000/  as  a  function  of  the  initial  value  plus  an  increment  due 
Af 

to  AC8. 

P 


Therefore 


dTooo/Aj 


PpP 

F 


1/2 


2/3 


C®  A 
P  f 


Thus,  we  can  write  the  expression 


(26) 


/V1* 

V  c  v, 


2/3 


1- 


2  A  Cp 


s 


3C 


i.  c. 


(27) 


where  l.o.  stands  for  the  Initial  condition  chosen  to  operate  about. 


p  p  ' 

Noting,  as  before,  that  J  F  ^  ^ 


f 


C  V 

p  f  i.c. 


we  can  write: 


T  T 

ooo  *  g 


‘f 


‘f 


2  AC 


1- 


3C 


i.  c. 


(28) 


This  is  a  more  convenient  form  since  T  is  easier  to  curve  fit  than  P_, 

g  F 


Letting  N  =  1  - 


%  g 
100 


we  find 

T  -  10,  000-62,  500N  +  125,  000N  lbs.  (29) 

g 

C  S=  .64  -  (.813  -  .000231  0  2)  R  -.  000069  0  2  (30) 

p  v  q  v 

-5  ,  2 

-2X10  0 


OR,  AC  S  -  (.  813-.  000231  0  2)  R  000069  0  2  -  2Xlo”50  2  (31) 
p  '  v  '  q  v  s 


8 

Once  was  computed,  it  was  possible  to  calculate  Tc  .  We  did 

^  s  s 

not  calculate  T  ,  but  the  value  1-T  .  This  is  equal  to  q/  s,  and  is 
c  c  q 

called  ”R  M.  AH  of  the  T  8  data  were  fitted  using  R  instead  of  T  S. 
q  c  g  q  c 

I?  is  desirable  because  at  T  =  1,0,  R  -  0,  and  is  more  accurately 
q  c  q  s 

calculated  on  the  analog  computer.  The  variations  of  with  T^  , 

/>  and  0  were  applied  through  a  circuit  with  a  0.  6  sec.  time  lag,  to 
V  s 

simulate  the  fan  dynamic  response. 


3.15.6  Nose  Fan  Force  and  Moment  Data 


The  same  technique  as  used  for  determination  of  main  fan 
Tqoo/a  was  applied  to  the  nose  fan  forces  and  moments. 

The  variation  of  nose  fan  lift  with  6  ,  thrust  reverser  door 

P 

position,  was  assumed  a  linear  function  of  6  ,  with  <5  =  0  correspond- 

P  P 

ing  to  maximum  lift.  The  nose  fan  maximum  normal  force  was  calcu¬ 
lated  as  a  function  of  airspeed,  N  ,  and  <5  ,  and  then  multiplied  by 

&  P 

<5p^?,  where  67°  <5^  was  the  point  of  maximum  negative  thrust.  Fig¬ 
ure  25  is  a  plot  of  per  cent  maximum  lift  versus  <5^. 

T  ,  which  is  the  thrust  of  the  pitch  fan  at  6  =0,  is 

°PF  P 

given  in  Figure  16 . 

The  lift  for  this  is: 

T  =  1500  +  .  1306V  +  .  0045V2  -  12100N  +  3400N2  (32) 
°PF 

and  using  the  previous  technique 

r  2ACS  1 


N„„  =  T  1 - 

PF  o  s 

max.  P„  3C 

F  P  < 


(Figure  30)  (33) 


AC  8  =  .000665<5  -4. 17X10~7<5  3  (6  in  deg.)  (34) 

P  PF  P  P  P 


N  =  N  max.  (1-.  0172  <5  )  +  3. 5w 
PF  PF  p 


The  pitch  fan  axial  force,  X 


XpF=  -  (4. 2-14N)  u  1- 


2  A  C, 
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The  moment  about  the  pitch  fan  hub,  Mn  ,  was  calculated  using 

°PF 


Figure  26  and  is  equal  to 


M  =  (1.766N  +  18.4N2)  Fi(V) 
°PF 

Where  F^  (V)  is  given  in  Figure  27. 


(37) 


The  effect  of  vectoring  on  the  6  values  is  shown  in  Figure  28. 
It  is  seen  that  the  nose  fan  6^  phases  out  to  zero  thrust  at  maximum 
vector  angle. 


Figure  29  shows  the  thrust  versus  N 
atlon  for  the  nose  fan  and  main  fans. 


g 


fer  single  engine  oper- 


Figures  31  and  32  are  the  computer  wiring  diagrams  for  the 
3°  transition  simulation. 


LONGITUDINAL  FORCE  COEFFICIENT 


NORMAL  FORCE  COEFFICIENT 


ANGLE  OF  ATTACK  -  I 


GROSS 

THRUST 


CG  AT  FS  246,  WL  112 
TAIL  OFF,  PITCH  FAN  OFF 


Figure  17  Estimated  Moment  Coefficient  About  C. 


CHANGE  IN  PITCHING  MOMENT  -AC 


RO  ANGLE  OF  ATTACK 
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PITCH  FAN  DOOR  POSITION  FROM  POINT  OF  MAXIMUM  LIFT  -  DEGREES 
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POSITIVE  LIFT 

=  25° 


0  =  20' 

v 


NEUTRAL  STICK 
POSITION 

(NEUTRAL  ELEVATOR) 


~j3  =  10° 
v 


FWD 


AFT 


■0  =  0° 
v 


PERCENT  OF  STICK  TRAVEL  FROM  MAXIMUM  FORWARD  POSITION 


Figure  28  Stick  Control  vs  Pitch  Fan  Door  Position 


*3 


2500  FT  HOT  DAY 


Figure  29  Estimated  Wing  and  Nose  Fan  Thrust  tor  Single  Engine  Operation 
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NOSE  FAN  POWER 


1500  LB  LIFT  SETTING 
HOT  DAY  -  2500  FT 


EST.  POWER  FOR  SINGLE 
ENGINE  OPERATION 


Figure  30  Nose  Fan  Power  Characteristics 
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4.0  LONGITUDINAL  THREE  DEGREE-OF-FREEDOM 
CONVERSION  SIMULATION 


4. 1  Discussion 

The  simulation  of  the  conversion  between  fan-powered  and  con¬ 
ventional  flight  was  utilized  for  several  purposes.  These  were: 

1.  To  provide  a  check  on  previously  calculated  maximum  fan- 
powered  flight  speed. 

2.  To  determine  requirements  of  tail  incidence  variation  during 
conversion  sequence. 

3.  To  develop  conversion  sequence  to  provide  smallest  transients 
at  conversion. 

4.  To  evaluate  handling  qualities  and  piloting  techniques  before, 
during,  and  after  the  conversion,  (below  power  off  stall  speed.) 

5.  To  evaluate  conversion  using  one  J-85  in  Jet  mode  and  one  in 
fan  mode. 

This  entire  study  was  conducted  using  /3V  =  50®,  <5f  =  45°,  ex¬ 
cept  that  later  runs  were  made  using  the  additional  15°  aileron  droop. 


The  method  of  approach  to  the  simulation  of  the  conversion 
aerodynamic  effects  deserves  mention. 

The  problem  of  generating  forces  and  moments  during  the  con¬ 
version  effectively  amounts  to  the  simulation  of  two  separate  airplanes. 
This  is  further  complicated  by  the  fact  that  during  fan  door  closure  the 
characteristics  are  for  neither  of  the  two  configurations. 

We  have  data  for  pure  fan-powered  flight,  and  also  data  for  con¬ 
ventional  flight  with  fan  closure  doors  open  and  closed.  We  also  have 
fan  run-up  and  run-down  data,  plotted  versus  time  from  diverter 
actuation. 
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The  problem  is  thus  presented  to  determine  how  to  determine  the 
force  and  moment  characteristic  at  conditions  intermediate  to  the  fan- 
powered  and  conventional  flight  regime. 

As  in  previous  work,  the  aerodynamic  data  was  separated  into 
conventional  data  plus  "fan  caused"  effects.  To  enable  this,  the  fan 
doors  open-fans  off  data  for  lift,  drag  and  moment  was  converted  to  slip¬ 
stream  notation  ,  and  subtracted  from  the  fans  on  data.  The  difference 
thus  obtained  was  plotted  as  a  function  of  Tcs  and  angle  of  attack.  These 
differences  were  then  the  net  effects  of  having  fan  power  applied  to  the 
aircraft.  As  fan  power  was  reduced,  the  differences  should  then  fade 
gradually  away  until  the  conventional  flight-fan  doors  open  data  remains. 

In  the  case  of  the  fan  to  conventional  conversion,  it  was  assumed  that 
during  the  fan  run-down  and  door  closing,  there  would  be  practically  no 
effect  on  fan  thrust  until  the  doors  were  almost  closed.  At  the  time  the 
doors  closed  fan  RPM  would  be  down  to  40  per  cent,  and  all  fan  effects 
would  be  switched  off. 

There  is  another  problem  concerning  the  conventional  flight  doors 
open  and  doors  closed  data.  There  was  no  intermediate  closure  door  posi¬ 
tion  data.  The  closure  doors  closed  data  was  subtracted  from  the  closure 
doors  open  data  to  obtain  increments  due  to  the  doors.  As  the  doors  went 
from  closed  to  open,  and  vice-versa,  the  increments  were  added  in  pro¬ 
portion  to  the  percent  of  door  motion  from  the  original  position.  These 
increments  are  relatively  small,  and  the  manner  of  phasing  them  in  and 
out  did  not  seem  to  be  critical.  Due  to  more  recent  fan  time  response 
data,  the  rate  of  door  closure  after  diverter  valve  actuation  was  slowed. 

The  first  use  of  tne  longitudinal  conversion  simulation  was  to 
check  the  previously  calculated  maximum  level  fan-powered  flight  speed 
attainable.  Previous  rough  calculations  had  indicated  the  aircraft  would 
trim  at  109  knots,  but  when  the  computer  was  turned  on,  the  vehicle 
slowed  down.  The  addition  of  the  J-85  ram  drag  aggravated  the  situation, 
and  for  these  reasons  the  new  wind-tunnel  data  for  45°  6  f  and  drooped 
ailerons  was  taken.  In  addition,  data  was  made  available  for  lower 
flight  speeds  (higher  Tc®). 

g 

At  about  this  time  it  was  determined  that  the  Cqo  values  obtained 
on  the  1/6  scale  fan-powered  model  were  higher  than  those  obtained  on 
the  full  scale  tests,  so  a  ACj)  8  correction  was  included,  derived  from 
comparison  of  full-scale  model  test  results  with  the  small-scale  data. 

The  drag  Increment  due  to  the  landing  gear  was  also  evaluated  so  that  the 
wheels-up  and  wheels-down  trim  speeds  could  be  determined. 
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Before  the  ACp  8  increment  was  included,  data  was  taken  from 
the  computer’ to  enable  cPetermination  of  trim  speeds.  For  this  purpose, 
the  J-85  RPM  was  set  at  100%,  and  the  lift  and  pitching  moment  were 
trimmed,  both  wheels-up  and  wheels-down,  for  various  velocities.  The 
resulting  net  horizontal  force  was  plotted  versus  velocity.  These  curves 
are  shown  in  Figure  33.  The  second  set  of  data  was  taken  at  100%  J-85 
RPM,  but  with  one  J-85  in  the  fan  mode  and  one  in  the  turbojet  mode,  for 
wheels-up  and  wheels-down.  A  large  increase  in  trim  forward  speed 
is  noted.  This  data  is  presented  in  Figure  34. 

The  third  set  of  data  have  the  ACn  8  correction  included,  and 
curves  are  presented  for  both  wheels-up  and  wheels-down  flight  in 
Figure  35. 

For  all  of  the  above  curves,  the  J-85  ram  drag  was  taken  as 
6. 25  Vy.  This  is  the  J-85  mass-flow  rate  at  2500  foot,  hot  day,  vel. 

«  200  fps,  multiplied  times  the  velocity.  This  value  is  conservative,  since 
some  of  this  drag  is  recovered  as  the  J-85  exhaust  is  directed  rearward 
by  the  vanes. 

To  determine  the  requirements  on  horizontal  tail  incidence  during 
conversion^  high-gain  autopilot  loop  was  connected  around  the  pitch  axis 
so  that  the  net  pitching  moment  would  be  zero  during  the  conversion. 

All  of  these  automatic  conversions  showed  that  an  essentially 
constant  tail  rate  was  required  as  the  butterfly  doors  cycled  through  their 
range. 


After  it  had  been  determined  that  a  constant  tail  rate  during  con¬ 
version  was  feasible,  pilot  controlled  conversions  using  various  conver¬ 
sion  sequences  were  performed. 

4. 2  Development  of  Conversion  Sequence 

After  it  had  been  determined  that  a  constant  tail  rate  was  accept¬ 
able  during  the  conversion,  there  remained  only  to  choose  a  sequence  of 
events  which  would  cause  the  least  conversion  transient. 

It  is  interesting  to  note  that  the  vehicle  gains  about  20  fps  in 
velocity  during  the  conversion  sequence.  This  enables  conversions  from 
fan  to  conventional  to  be  started  below  the  stall  speed. 

The  following  conditions  were  obtained  prior  to  conversion  for  eval¬ 
uation  of  the  conversion  sequencing:- 


I 
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Prior  to  Conversion  to  Fan  Mode 

1.  Flaps  down  full  (45°),  aileron 
drooped. 

2.  Landing  gear  down. 

3.  fi  v  =  50° 

4.  Nose  fan  inlet  aud  exit  open. 

5.  Thrust  spoilers  out. 

6.  100%  Ng. 

7.  Vtrim  =  203  fps  EAS 


Prior  to  Conversion  to  Conventional 
Mode 

1.  Flaps  down  full,  ailerons 
drooped. 

2.  Landing  gear  up. 

3.  p  v  =  50°. 

4.  Nose  fan  inlet  and  exit  open. 

5.  Thrust  spoilers  retracted. 

6.  Ng  =  100%. 

7.  vtrim  =  170  EAS 


The  objective  was  to  find  a  conversion  sequence  which  would 
cause  the  least  transient  and  result  in  the  minimum  amount  of  pilot  effort. 


During  the  fan  to  conventional  conversion,  the  aircraft  must 
change  attitude  from  a  slightly  negative  angle  of  attack  up  to  about  a  ten 
degree  positive  angle  of  attack.  In  the  conventional  to  fan  conversion, 
the  rotation  will  be  less  because  the  vehicle  is  at  a  smaller  initial  angle 
of  attack  due  to  higher  speed.  The  sequences  picked  result  in  a  practically 
hands-off  conversion,  and  changes  from  the  desired  flight  conditions  only 
require  small  pilot  corrections. 


The  desired  conversion  sequences  were  determined  to  be  as 

follows: 


Function  Conv.  to  Fan 

1.  Diverter  valve  delay  time.  .  8  Sec. 

2.  Diverter  Valve  Actuation  .  4  Sec. 

time. 

3.  Fan  door  delay  time.  0 

4.  Horiz.  Stab.  Delay  Time.  .  5  Sec. 


Fan  to  Conv. 
0 

.4  Sec. 


.4  Sec. 
0 


(Continued) 


0 


5. 

Fan  door  actuation  time. 

1. 5  Sec. 

3.0  Sec. 

6. 

Horiz.  Stab,  actuation  rate. 

7. 5°/Sec. 

7. 5°/Sec. 

7. 

Final  Horiz.  Stab,  position. 

o 

00 

+ 

-5° 

8. 

Horiz.  Stab,  travel 

12.6° 

16.6° 

In  general  the  conversion  was  easily  accomplished  as  long  as  the 
horizontal  tail  moved  in  the  correct  direction,  but  as  the  tail  sequence 
departed  from  the  optimum,  more  pilot  control  was  required. 

The  conversion  from  fan  to  conventional  flight  was  also  evaluated 
at  an  initial  velocity  below  the  power-off  stall  speed.  This  conversion  is 
practical  because  during  the  3  seconds  of  the  conversion  sequence,  the 
unspoiled  J-85  thrust  results  in  an  average  acceleration  of  about  10ft/sec2, 
for  a  velocity  gain  of  30  fps  or  17. 75  knots. 

Care  must  be  taken,  however,  because  too  large  an  angle  of 
attack  can  result  in  sufficient  drag  to  result  in  no  velocity  increase. 

Figure  36  is  a  recording  of  several  variables  during  a  nominal 
fan-conv.  conversion  (both  hands-off  and  pilot-controlled). 

Figure  37  is  a  recording  of  a  nominal  conv. -fan  conversion. 

Figure  38  and  39  -  Single-engine-piloted  -  both  ways. 
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Figure  35  Net  Horizontal  Force  v/lth  Trim  Lilt  and  Moment  vs  Velocity  lot  Both 
Engines  in  Fan  Mode 


4.3  Analysis 


Equations  of  Motion 

The  conversion  maneuver,  being  restricted  to  small  angles  of 
attack,  was  simulated  using  a  set  of  longitudinal  stability  axis  equations 
of  motion.  The  variables  in  the  stability  axis  equations  are  total  velocity, 
flight  path  angle,  and  pitch  angle.  Figure  40  shows  the  stability  axis  coor¬ 
dinate  system  and  the  three  coordinates. 

It  is  seen  that: 

y  =  q  -  a  (*) 


X  and  Z  are  the  body  axes,  rotated  through  the  angle  a  from  the 
velocity  vector. 

The  lift  is  perpendicular  to  the  velocity  vector,  and  the  drag  is 
parallel  to  the  velocity  vector. 


We  can  thus  write  the  equations  of  motion: 


g  sinT 


D_ 

m 


a 


gcosy  | 


;  sm 

y 

Where: 


(2) 


(3) 

(4) 


D  * 


s  s 

CD  q  A 


f 


D  „  =  Nose  fan  drag 
NF 

Tj  =  Turbojet  thrust 
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s  s 

L  =  C  q  A, 
L  I 


L  =  horiz.  tail  lift  =  .0486 

l 


i.  -4°  +  4t 


qs, 


s  s 


Zm  =  C  q  A,D,  -  .02 
m  ii 


.  .  o  6  6 

lt'4  +  ~T 


qS  c 
w 


*  42. 6  ft. 


Df  =  5.22  ft 


st  = 


50.6  ft 


S  =  260  ft 
w 

P  =  .  00205  slugs/cu.  ft 


m  --  9200/g  slugs 
,  2 

g  =  32.  2  ft/ sec 
c  =  9. 4  ft 


€  =  downwash  at  tail  at  zero  angle  of  attack,  degrees. 


T  ,  Wing  fan  Lhruot  loading,  psf. 

ooo/A^. 


q  =  q  +  i 


lino/ A, 


ix  =  tail  incidence,  deg. 
t 

P  2 
Q  “  2  VT 


s 


Elevator  Pos.  deg. 
T 

ooo/A. 


q  4  T 


ono/A, 


66Vm0 

T 
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Figure  41  Eltect  ot  Wing  Closure  Doors  -  Power  Oil 
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Figure  45  Thrust  vs  Time  irom  Diverter  Actuation  tor  Nose  Fan  and  Main  Fan  flunup 

V  -  100  KNOTS 


P 

i 

a 

2 

Q 

; z 

< 

w 

8 

Z 


TIME  -  SECONDS 

Figure  46  Variation  of  Nose  Fan  Ram  Drag  Following  Diverter  Valve  Operation 
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VI 

Lfh.ijith 

§ 

s 


R  =  1  -  T 
q  c 


s 

C  =  Slipstream  Jilt  coefficient 

1  (includes  tail  lift  for  i  =  4°) 

v 

g 

Cp  =  Slipstream  drag  coefficient 


Fig.  41  and  Fig.  42 

Fig.  41,  Fig.  42  and 
Fig.  43 


Slipstream  Pitching  moment  coefficient  Fig.  44  and  Fig.  41 
(includes  tail  moment  for  i  =4°) 

v> 


The  set-up  of  this  data  required  the  simulation  to  cover  the 

g 

T  range  of  about  .  9  to  zero,  wing  fan  doors  open  and  closed, 
c 

The  curves  of  Figure  41  depict  the  \  riation  in  the  longitudinal 
aerodynamics  due  to  the  position  of  the  wing  fan  closure  doors,  in  the 
power  off,  or  Tc  =0,  configuration. 

The  curves  of  Figure  42  and  44  show  the  variation  of  the  longitudinal 
aerodynamics  with  the  wing  fan  closure  doors  open  and  at  various  T  s 
values  applicable  to  the  conversion  regime. 


In  addition,  the  curves  of  figures  45,  46,  47,  and  48  were  used  to 
simulate  the  effect  of  fan  speed-up  and  run-down  during  the  conversion 
maneuver. 

g 

The  data  was  simulated  by  subtracting  the  doors  open,  T  =0,  data 

of  Figures  42  and  44.  This  remainder  was  considered  to  be  the  fan-caused 

effects,  and  the  effect  of  fan  run-down  was  applied  to  these  increments, 
s  s  s 

These  AC^  ,  and  AC^  ,  and  AC^  are  shown  in  Figures 49,  50  and  51, 
plotted  at  constant  values  and  zero  angle  of  attack. 

s  s  s 

(Example)  The  component  of  C  from  T  =  0  reflected  to  the  Tc  =.954 
curve  on  Figure42is  equal  tothe  C I  T =  0  value  multiplied  by 
s 

(1-T  ).  These  differences,  then,  are  plotted  in  Figures  49,  50  and  51.  As 
the  fan  thrust  decays  during  conversion,  T„  changes  and  the  increments 
due  to  fan  operation  automatically  phase  out.  When  the  doors  close,  the 
remaining  small  increment  is  simply  turned  off. 


During  the  door  closing  (or  opening)  period,  the  increments  ob¬ 
tained  from  Figure  41  are  phased  in  or  out  proportional  to  door  position. 

The  3°  conversion  longitudinal  computer  diagram  is  shown 
in  Figure  51A. 
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5.0  SIX-DEGREE-OF-FREEDOM  PERTURBATION 
CONVENTIONAL  FLIGHT  SIMULATION 


5. 1  Discussion 


The  simulated  handling  qualities  of  the  XV- 5A  aircraft  were 
investigated  for  a  range  of  flight  conditions  and  aircraft  configurations 
considered  representative  of  the  conventional  flight  regime.  The  aero¬ 
dynamic  data  used  represent  the  present  best  estimate  of  the  aircraft 
aerodynamic  characteristics  based  on  wind  tunnel  data  and  standard  es¬ 
timating  procedures  where  data  were  lacking. 

Handling  qualities  were  e  valuated  using  the  Cooper  rating 
system  and  the  results  are  presented  in  Tables  I  and  II  following.  Char¬ 
acteristics  were  investigated  for  gross  weights  of  7500  and  12,  500  pounds; 
center  of  gravity  at  fuselage  stations  240  and  246  inches;  altitudes  of  sea 
level,  20,  000  and  32,  000  feet  and  a  speed  range  from  85  to  500  knots 
(M  =  .  13  to  .  84).  Table  I  presents  a  complete  list  of  conditions  studied 
and  Table  II  gives  the  evaluation  of  each  condition  for  a  series  of  maneu¬ 
vers  in  terms  of  Cooper  pilot  ratings. 

An  overall  evaluation  of  the  XV-5A  flying  qualities  in  the 
conventional  flight  mode  would  be  Cooper  Rating  3,  "satisfactory, "  but 
with  some  mildly  unpleasant  characteristics.  The  "mildly  unpleasant" 
refers  to  weakening  dihedral  effect  during  conditions  of  zero  or  slight 
negative  angle,  such  as  conditions  2  and  4,  flaps  down,  180  knots  and 
flaps  up  500  knots  respectively.  The  problem  appears  when  the  aircraft 
develops  a  slight  left  bank  while  developing  a  right  yawing  velocity,  due 
to  a  right  rudder  input  in  the  absence  of  co-ordinating  aileron  during  the 
maneuver  "rudder-S-turns". 


This  problem  is  in  some  degree  a  result  of  adverse  roll  due 
to  rudder  Cg^  which  was  verified  by  removing  C£^  and  thereby  achiev- 

satisfactory  dihedral  effect.  Cg^  for  angle  of  attack  zero  is  -.  00127/deg. 

and  C/,  is  .00025/deg.  The  pilot  rates  this  specific  maneuver  between 
6 

r 

5  and  6  "unacceptable  for  normal  operation"  and  the  condition  rates  4. 


in 


This  lack  of  dihedral  effect  at  these  conditions  is  not  apparent  during  any 
of  the  coordinated  maneuvers. 


Low  speed  flaps  down  or  condition  1  was  studied  extensively. 
The  overall  rating  for  condition  1  was  3.  Considering  the  flight  speed 
of  85  knots  or  15%  over  stall  speed,  the  aircraft  "flies  well".  To  further 
investigate  this  condition  the  longitudinal  stability  was  increased  and  de¬ 
creased  by  varying  Cm  and  Cm  ,  It  was  found  that  even  large  changes 

Ot  Q 

in  C  or  C  (50%)  had  little  effect  on  the  pilot's  ability  to  fly  the 
ma  mq 

simulator.  Reducing  C  or  C  to  zero  made  less  than  1  change  in 

iii/u  m 

or  q 

Cooper  rating.  However,  reducing  both  £m  and  Cm  had  a  much  more 

deleterious  effect  on  handling  qualities  and  Cooper  rating  jumped  to  6. 

This  study  indicates  that  any  loss  in  static  longitudinal  stability  associated 
with  flight  at  low  speeds  and  high  angle  of  attack  would  not  be  a  threat  to 
flight  safety. 


At  the  high  jpeed  condition  some  studies  were  conducted  to 
establish  the  aircraft's  sensitivity  to  certain  parameters  which  might 
change  as  a  means  of  improving  handling  qualities  or  due  to  possible  errors 
in  the  aerodynamic  data  used.  Increases  in  C ^  of  50  and  100%  greatly 


improves  the  handling  characteristics  particularly  "wing  dropping"  tenden¬ 
cies.  Increases  of  this  magnitude  in  C  do  not  greatly  change  the  han- 

n/j 

dling  qualities.  Very  large  changes  in  C  ,  yawing  moment  due  to  ai- 

n<5 

a 

leron,  (as  much  as  1000%)  in  either  direction  positive  or  negative  had 
little  or  no  effect  on  general  handling  qualities. 


Changes  in  altitude,  gross  weight,  and  center  of  gravity  did  not 
greatly  affect  handling  qualities.  For  example,  changes  to  high  altitude  or 
heavier  weight  were  evaluated  as  "same  as  before  or  maybe  a  little  less". 
That  is  to  say,  if  a  maneuver  got  a  Cooper  rating  of  3  at  a  given  condition, 
it  would  not  vary  more  than  one-half  for  other  altitudes,  gross  weights,  or 
centers  of  gravity.  For  this  reason  and  due  to  time  considerations,  some 
conditions  which  were  originally  scheduled  for  checks  were  eliminated. 


5.2 


Summation 


1.  Reduced  static  longitudinal  stability  during  low  speed 
flaps  down  high  angle  of  attack  flight  does  not  present 
a  handling  qualities  problem. 

2.  A  50%  change  in  directional  stability  does  not  greatly 
alter  the  overall  handling  qualities. 

3.  A  50%  increase  in  lateral  stability  would  improve  fly¬ 
ing  qualities  considerably. 

4.  Overall  flying  qualities  of  the  XV- 5 A,  as  is.  are 
satisfactory. 

5.  Changes  in  aircraft  characteristics  which  result  from 
changes  in  altitude,  gross  weight,  and  center  of  gravity 
do  not  greatly  alter  handling  qualities. 

*  For  the  purpose  of  this  repo  *,  a  change  in  handling 
qualities  which  changes  Cooper  rating  by  more  than 
one  is  large,  whereas  less  than  one-half  is  small. 
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5.3 


Analysis 


The  following  are  the  standard  small  -  perturbation  lateral- 
directional  and  longitudinal  three-degree-of-freedom  uncoupled  equations 
of  motion  used  throughout  the  industry. 

5.3.1  Longitudinal  Perturbation  Equations 

The  longitudinal  perturbation  equations  are: 


u=X  u+X*w+X  w  +  x  0  +  X  6  -g(cosrn)0 
u  w  w  q  6e  e  ° 


w=Zu+Z*w+Z  w  +  (U  +  Z  )0  +  Z  6  -  g(sin  7  Q) 
u  w  w  o  q'  dee 


V  =  Mu+M.w  +  M  w+Mfi+M,  6 
u  w  w  q  6e  e 


Changing  all  angles  to  degrees,  and  letting  w  =  U<a  /57. 3 

u  =  A  u  +  ..a+Aa  +  A*10  +  A  6  -  JL  (cosyo)0 
u  a  a  e  a0  « 


a  =  B  u  +  B.  a  +Ba  +  B.  0  +  B.  6  -  JL  (sin  V  ) 
u  a  a  q  (5e  e  u  o' 


(1) 

(2) 

(3) 

(4) 

(5) 


0  =  C  U  +  cat  +  C  a  +  C*  0  +  C,  6 
u  a  oi  e  fie  e 


(6) 


Where 


A  =  X 
u  u 

Aa  =  XwV57*3 
Ao=XwV57-3 

A*  =  X  /57.3 

0  q 

A.  =  X,  /57.3 
o  6 


B  =  Z  (57.3/U  ) 
u  u  o 

B.  =  Z* 


a  w 

B  =  Z 
a  w 


B  •  =  1  +  Z  /U 
e  q  o 


e 


e 


C  =  57.3  M 
u  u 


C.  =  U  M* 
ot  o  w 


U  M 
o  w 


C-  =  M 

0  q 


e 


lit 


o  t? 


5.3.2  Lateral  Perturbation  Equations 


>  =  V  +  V  +  V  +  Yaa*  a  +  YaA  ‘  V* 


-*♦+ 


+  7jt(sin1'o)* 
o 


*  =  V  +  V  +  V  +  Laa4a  +  Larar  +—  * 


*  = 


p  p  r  oa  a  or  r 


(7) 

I 

xz 

I 

•  • 

(8) 

X 

I 

0  » 

xz 

(9) 

5. 3. 3  Determination  of  and  C 

u  u  6  u 

e 


L-CL  <a-c*oL)  +  CL  6e  +  2U<Cl“  +  CL  *> 

q 


a 


6e 


a 


(10) 


CL  =  f  +  fCT  .  i  +  fC 
D  1  2  L  fie  I  3 


ijae| 


Cm  =  Cm  +  Cm  “  +  Cm  a»  +  2U(C  „  +  C  9  ) 


m  m  m  me 

o  a  6e 


m .  m 
a 


(11) 

(12) 


where 

f. 


f  (C  ,  M)  is  the  untrimmed  drag  polar 
1  Li 


f  =  i  (M)  is  the  trim  drag  coefficient 

u  4 


D<5e 


f  -  f  (M)  is  the  trim  drag  coefficient 

O  U 


'D6e 
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5.3.4  Dimensional  Stability  Derivatives  in  Terms  of 
Non-Dimensional  Stability  Derivatives 


Longitudinal 


x  -  -£su  (-c  -  c  ) 

U  m  D  °u 


Z»  ’  ^'-CL  -  CL  1 
m  u 


M  «  ^?(C  +  C  ) 

u  ly  m  mu 


+  -  (4^)  ^08(0  +  »  ) 
m  '  8u  T  o 


7T  8T 


PSU 

X  =  Sm  (CT  -  ) 

w  '  L 


*.  •  ^  (‘CL  '  Ct) 


(Cm  > 

y  m  a 


x  =-^-  (-C  ) 

6  2m  D 
e  fle 


Z .  ,  -H*L  <-C  ) 

w  4m  L  • 


M  .  (C  ) 

w  4i  m  . 

y  “ 


z  ,  £S£(-cf  ) 

q  4m  Lq 


m  .  -£M2  (c  . ) 

q  41  m’ 

y  * 


\  <-CL  > 

fle  2"  \ 


n,  ^h1c(C  , 

21  m . 


Lateral 


PSU 

Y  =  -  C 

P  2m  y. 


PSU2b  „ 

L  - -  co 

’  a,  ‘f 


PSU2b 

N*  ■  ~  % 


PSb 

Y.  - ~~  C 
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«  41  n. 
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PSUb 
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PSUb 
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PSU  „ 

\  2m  V; 


.  PSU2b 
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Definitions 


T  Thrust,  lbs 

ot  ^  Angle  between  thrust  axis  and  body  x-axis 

aQ  Initial  angle  of  attack 

z^,  Thrust  moment  arm,  ft 


5.3,5  Conversion  of  Inertia  Data  From  Body  Axes 
to  Stability  Axes 


If  Ixb*  Izb»  and  Ixzb  are  the  moments  and  product  of  inertia 

expressed  in  body  axes,  and  <r  is  the  angle  between  the  body  x-axis 

and  the  stability  x-axis,  then  the  moments  of  inertia  and  product  of 

inertia  in  the  stability  axis  system  (1,1,  and  I  )  are: 

X  z  xz 

'x  =  ('zb  +  W/2  -  (Izb  -  W  (COS  2  a)/2  •  Ixzb(,in  2  *>  <18> 


lz  =  (Izb  +  W/2  +  (I-K  -  W  «*»  2* >/2  + 


zb  xb 


xzb 


(19) 


!XZ  =  Ixzb<C0S2ff)  -  ('zb  '  W  (Sin  2  > 


(20) 


NOTE:  If  conventional  body  axes  are  used  a  =  a 


■  *  ►  STAI 


STABILITY  X-AXIS 


Equations  (1),  (2),  (3),  (7),  (8),  and  (9)  were  modified  to  include  coup¬ 
ling  terms  from  the  lateral  to  the  longitudinal  modes.  In  addition,  the 
gravity  components  were  computed  using  the  Yaw- Roll- Pitch  direction 
cosine  relationships.  The  gravity  perturbation  values  were  then  com¬ 
puted  by  subtracting  the  gravity  values,  obtaining  at  the  initial  flight 
altitude. 


The  modified  equation  of  motion  became  as  follows: 

u=Xu+Xw+X,  6e  +  (g  -  g  )  +  rv-  qWo  (21) 

u  w  5e  x  xo 


w  =  Z  u  +  Z.w  +  Z  w  +  Zq  +  Z,  +6e  +  (g  -  g  )  +  qU  -  pv  (22) 
u  w  w  q  oe  z  zo  o 


*  =  V  +  Y« aaa  +  V  +  V  '  r  +  (sy  '  gyo>  +  p“  o  <26> 


Where: 


and 


=  g l  =  -  g  Sin  0  Cos  <p 

Z  0  0 

o 

=  grng  =  g  Sin 


=  gn^  =  g  Cos  Cos  <£o 


0  =  r 

0  =  P 

6  =  q 


12o 


/ 


=  initial  x  axis  velocity  component 
Wq  =  initial  z  axis  velocity  component 
H^x  axis  component  of  power  plant  angular  momentum 

X 

=  z  axis  component  of  power  plant  angular  momentum 

The  computer  diagrams  used  to  mechanize  these  equations  are  shown 
in  Figures  52,  53  and  54. 
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6.0  QAS  GENERATOR  CONTROL  ANALYSIS 


6. 1  Discussion 


The  possibility  of  using  J-85  RPM  alone  for  hovering  thrust 
control  was  investigated  concurrently  with  the  hovering  simulation.  The 
fan  response  was  approximated  by  an  equivalent  first  order  lag  of  r  =  .  6 
seconds  at  the  fan  speed  necessary  to  enable  hovering  at  2500  feet  on  a 
hot  day.  The  J-85  response  was  simulated  using  a  small-perturbation 
analysis  around  98%  J-85  RPM,  and  including  the  limiting  effects  of  th^ 
fuel-control  system.  The  computer  wiring  for  the  J-85/  Fan  combination 
is  shown  in  Figures  56  and  57  following. 

Let  us  assume  for  the  moment  that  the  combined  engine-fan 
response  can  be  approximated  by  the  transfer  function 
Thrust  ~  1 

Thrust  Command  (1  +  .  6S)(1  +  t^S  )  the  .  6  second  time  con¬ 
stant  being  the  fan  equivalent  fan  lag  and  the  engine  being  represented 
by  a  first  order  lag  of  r e  .  This  is  not  a  true  representation  of  the 
engine,  but  it  will  suffice  for  this  discussion. 

If  we  can  introduce  a  transfer-function  between  the  throttle 
and  engine  of  1  *  TjS  ,  it  would  be  possible  to  make  t  ^  =  .  6  second, 

1  +t2S 

and  thus  cancel  the  engine  .  6  second  lag. 

E0  (1+KtS) 

This  transfer  function  may  be  written  also  as  - —  =  k. - - - , 

eIN  1  1  +  t  s 

where  r  is  fixed  at  an  arbitrary  value  and  Kg  is  called  the  "magnification 
factor".  The  response  of  this  transfer  function  to  a  step  input  shown  in 
Figure  57,  assuming  =  1.0  and  K2  =  4. 

An  example  of  an  electrical  network  to  realize  the  transfer 
function  is  also  shown  in  Figure  57. 
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This  network  has  the  transfer  function 


where  we  can  let  K. 


K  .  h^L 

«2 


+  R2 


in 

RoR.  C  . 
2  1  ,  and 

RlR2 


H2  (1  +  RjCS) 
*1*2 


/l  +  R.RnCS  < 
(  R1 +  R2  / 


A  passive  hydraulic  unit  could  be  mechanized  to  give  the  same 
transfer  function. 

By  applying  this  transfer  function  between  the  throttle  and  fuel 
control,  we  obtain  the  complete  system  transfer  function 
Thrust  _  Kx  (1  +  K2  t  S) 

Thn,St  pilot  command  =  V-WXl+V*  (1+tS).  It  is  seen  that  we 

can  choose  K2  and  r  so  that  the  1  +  .  6S  term  is  cancelled,  and  the  result¬ 
ing  transfer  function  has  the  engine  lag  plus  a  new  lag  which  has  a  shorter 
time  constant  than  the  fan  lag.  In  this  manner  the  system  response  can 
be  improved. 


However,  there  is  a  limit  on  how  large  K2  can  be,  and  this 
therefore  limits  the  value  of  t  ,  Due  to  the  fuel-flow  limiting  effects 
of  the  acceleration  cam  in  the  J-c  fuel  control,  increasing  K2  to  a  large 
value  does  not  necessarily  result  in  improved  system  response. 

Assuming  98%  gas  generator  RPM  is  required  for  hovering  flight 
at  2500  feet  hot  day  conditions,  two  separate  conditions  of  fuel-flow  limit¬ 
ing  were  investigated.  A  set  of  data  was  taken  with  a  100%  N  g  throttle 
limit,  using  a  step  1%  A  Ng  command  starting  at  98%  N  .  Another  set  of 
data  was  taken  under  the  same  conditions  but  using  a  102%  N  throttle 
command  limit.  The  parameter  measured  was  the  time  for  uie  thrust  to 
rise  to  90%  of  the  final  thrust  change,  with  the  restriction  that  the  over¬ 
shoot  should  be  less  than  20%  of  the  thrust  change.  Figure  58  is  a  plot  of 
settling  time  VS  magnification  factor  and  r  for  100%  N  limit, 

^command 

and  Figure  59  is  the  same  plot  for  a  102%  N  limit. 

^command 
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The  improvement  in  response  tapers  off  as  tne  magnification 
factor  goes  up,  because  an  increased  input  to  the  fuel  control  does  not 
result  in  an  increase  in  fuel  flow  once  the  throttle  oommand  limit  or 
acceleration  cam  unit  is  reached.  Thus,  beyond  a  reasonable  value,  an 
Increase  in  the  magnification  factor  does  little  good.  If,  however,  the 
gas  generator  is  able  to  operate  well  below  the  acceleration  cam  limit, 
such  as  at  8.L.  standard  day,  Increased  magnification  can  be  utilized  to 
advantage.  It  is  desirable  to  limit  the  engine  speed  by  limiting  the  throt¬ 
tle  input,  because  excessive  use  of  the  acceleration  cam  limit  can  cause 
a  decrease  in  engine  life  due  to  higher  average  turbine  temperature. 

As  a  result  of  the  computer  runs  used  to  develop  Figures  58 
and  59 ,  a  r  of .  15  second  and  a  magnification  factor  of  5  have  been 
chosen  as  optimum.  The  time  to  90%  of  final  value  is  reduced  from 
1. 4  second  to  .  6  and  .  5  seconds  for  throttle  command  limits  of  100% 
and  102%  respectively. 

Representative  computer  records  are  shown  in  Figures  60,  61 
and  62  for  the  bare  engine  and  fan,  and  for  the  recommended  jazzer  values 
and  100%  and  102%  throttle  input  limits. 

In  addition  to  the  computer  runs,  which  include  the  non-linear 
effects  of  the  fuel  control  and  throttle  input  limit,  the  linearized  block 
diagram  furnished  by  General  Electric  was  solved  for  the  small-pertur¬ 
bation  linear  J-85  response. 


The  J-85  block  diagram  and  small  perturbation  analysis  ar« 
given  in  Section  6. 2,  with  pertinent  engine  data. 


The  J-85  closed-loop  response  for  small  inputs  at  100%  RPM 
was  thus  derived  as 


55. 1  (S/100  +  l)/s2  1. 8S/ 

— -  +  '29  +  1 

_ \292 _ 

'jl  +  V-iL  ♦  lie  +  +L84s 

^180  /  \  99. 72  99,7  /  \46. 32  46,3 

/  S2  +  .  986S  +  A 

\  12. 242  12.24  / 


This  hand  analysis  was  used  to  check  the  computer  mechanization  of 
the  85* 
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Figure  60  Computer  Output  Records  lor  Bare  EngiaePlus  Fan  Response  lot  1%  A  Ng  Input 
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Figure  61  Computer  Output  Records  tor  100 %  Figure  62  Computer  Output  Records  for  102 % 

Ng  Limit  and  1%  A  Ng  Command  jVg  Limit  and  1 %  A  Ng  Command 


6.2 


Analysis 


The  data  used  for  the  small-disturbance  hand-calculation 
of  J-85  response  are  given  in  Sections  6. 2. 1  and  6. 2. 2.  Figures  63 , 
64,  and  65  are  plots  of  the  data  of  Section  6. 2. 1.  Figure  66  is  a  plot 
of  the  95°  F,  2500  foot  altitude  acceleration  cam  limits  on  ^f»  as  well 

P T 

as  the  steady  state  operation  line  from  Section  6.2. 1. 


S3 


Figure  67  is  a  copy  of  the  original  block  diagram  of  the  engine 
and  fuel  control  system  as  furnished  by  General  Electric.  The  following 
Figures  68  through  73  show  the  successive  reductions  of  this  block  dia¬ 
gram  to  a  single  block,  giving  the  linear  closed-loop  J-85  response  to 
throttle  inputs. 

6.2. 1  Wj  and  Pgg  Values 
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UPPER  LIMIT  ACCEL.  SCHED.  (95 


Figure  66  Operational  Line  and  Schedule  lor  2500'  Hot  Day 
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(bTHAjS-'- 

Figure  67  Block  Diagram  of  J-8S  Engine  and  Fuel  Control  System 


Figure  68  Reduced  Block  Diagram  ot  JSS  Engine  and  Fuei  Control  System 


Figure  69  Further  Reduced  Block  Diagram  ot  J-85  Engine  and  Fuel  Control  System 
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Figure  70  Further  Reduced  Block  Diagram  o t  J-8S  Engine  and  Fuel  Control  Syste 
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Figure  71  Further  Reduced  Block  Diagram  oi  JSS  Engine  and  Fuel  Control  System 
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Figure  72  Final  Engine  Closed-Loop  Small  Disturbance  Transfer  Function 
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Figure  73  Final  Engine  Closed-Loop  Small  Disturbance  Transfer  Function  (Factored) 


6.2.2 


Engine  Partials  and  Control  Data 


Engine  Data 
N-  Constant 


A. 

Standard  Dav  S.  L.  S. 

B.  2500'  Hot  Day 

TT2  *  »3' 

PT2  -  13.  5 
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»0%N 

95%N 

100%N 

90%N 

95%N 
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f  /  awf 

.37 
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Available  via  J-85  deck  at  Evendale 


Fvt  =  Net  Turbojet  Thrust 


100%N  •  16,  500  RPM 


J-  16.5  lb  -  ft 


Control  Data 


Constants 

$  s.s. 

fcf 

.9 

Damp.  Ratio 
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ft 

29  rad/sec 
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< 
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(Wj/PgJmin 

= 
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(Wf)min 

= 
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TP 
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P3  Serv  o  Lag 
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Reduction  of  Block  A 
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_ 55. 1  a  +  ,01S)  (S2/292  +  1.88/29  +1) _ 
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(S/180  +  1)  (S2/99.72+  1.76S/99.7+  1)  (S2/46.32+  1.84S/46.3+  1) 
(S2/12.  242+  .  986S/12.  24+  1) 


7 . 0  LIFT-  FAN  ROLL- YAW  CONTROL  COUPLING  ANALYSIS 

7*1  Discussion 

Coupling  from  roll  to  yaw  and  yaw  to  roll  is  induced  for  con¬ 
trol  inputs  at  other  than  zero  collective  vector  angle.  This  effect  is  most 
serious  in  the  yaw  mode,  where  a  yaw  control  imput  can  result  in  a  larger 
rolling  moment  than  yawing  moment.  Coupled  with  the  fact  that  the  roll 
moment  of  inertia  is  about  1/5  the  yaw  moment  of  inertia,  it  is  seen  that 
yaw  control  imputs  can  cause  large  roll  accelerations. 

This  effect  has  been  compensated  for  in  the  mechanical  mixer 
located  between  the  control  stick  and  the  exit  louver  inputs.  Yaw  due  to 
roll  input  is  not  compensated,  because  the  effect  is  less,  due  to  the  higher 
yaw  inertia,  and  furthermore,  the  effect  is  in  the  direction  of  the  roll,  so 
that  a  right  roll  induces  a  right  yaw.  The  roll-due-to  yaw,  on  the  other 
hand,  is  in  the  adverse  direction. 

The  compensation  of  roll  due  to  yaw  consists  of  cross- coupling 
the  yaw  and  roll  signals,  while  gradually  phasing  the  yaw  louver  inputs  to 
zero  at  45°  (3  .  At  this  vector  setting,  all  yaw  control  moments  are  caused 
by  rudder  deflections  alone,  as  far  as  pilot  inputs  are  concerned. 

As  far  as  the  SA  system  is  concerned,  no  mixing  of  yaw  and  roll 
signals  has  been  decided  upon,  since  the  effects  of  this  cross-coupling  (in 
the  SA  system)  on  flying  qualities  has  not  as  yet  been  determined.  If  the 
system  behavior,  as  determined  below  is  deemed  objectionable  to  the  pilot, 
then,  either  an  electrical  mixer  and/or  phasing  out  of  yaw  will  be  necessary. 

The  effect  of  vectoring  on  the  roll  and  yaw  moments  caused  by 
differential  stagger  and  differential  vector  is  shown  in  Figure  74.  Essen¬ 
tially,  as  the  collective  vector  angle  increases,  the  effective  axes  of  moments 
for  roll  and  yaw  inputs  rotate  with  the  vector  angle,  and  approximately 
proportionally.  Fortunately,  this  rotation  is  in  such  a  direction  that  it  results 
in  a  stable  system.  Rotation  in  the  opposite  direction  can  result  in  a  diver¬ 
gent  system. 

The  rate  gyros  are  fixed  in  the  vehicle,  and  thus  measure  the 
angular  rates  in  an  axis  system  displaced  from  the  axis  of  application  of 
corrective  signals. 


As  is  shown  in  Section  7. 2  following,  a  simple  inertial  two- 
degree-of-freedom  system  was  investigated  which  allowed  only  yaw  and 
roll  angles  to  occur.  No  translations  were  included.  This  two-degree- 
of- freedom  representation  results  in  a  4th  order  characteristic  equation, 
which  has  as  a  solution,  two  sets  of  second  order  roots.  One  pair  of  roots 
would  be  the  "roll"  mode,  and  the  other  pair  the  "yaw"  mode. 

At  35°  /3  the  characteristics  of  the  response  have  been  deter¬ 
mined  for  roll  gains  of  1  and  10,  using  a  compensation  ratio  setting  of  .  1 
for  both  axis.  The  yaw  gain,  K  ,  was  varied  over  the  range  of  1  to  20. 

y 

Figure  75  is  a  plot  of  the  damping  ratios  and  natural  frequencies 
for  =  2,  and  variable  K^. 

Figure  76  is  the  same  plot  with  roll  gain,  K  =  20. 

K 

It  can  be  seen  that  the  variation  of  the  yaw  gain  results  in  a 
variation  of  the  characteristics  of  both  the  roll  and  yaw  modes,  since  both 
natural  frequencies  vary. 

This  is  not  desirable,  but  may  be  tolerable.  This  must  be  deter¬ 
mined  from  pilot  evaluation  tests. 

It  would  be  most  desirable  to  have  the  roll  gain  affect  only  the 
roll  axis,  and  the  yaw  gain  the  yaw  axis.  In  this  manner,  the  handling 
qualities  for  a  given  se<  of  roll  and  yaw  gains  would  be  unique  and  determin¬ 
able.  If  this  coupling  is  present,  the  handling  qualities  for  a  given  set  of 
gains  would  not  necessarily  be  unique,  and  it  is  very  likely  that  essentially 
the  same  characteristics  could  be  arrived  at  with  several  different  gain 
combinations. 

With  this  in  mind,  and  realizing  that  the  axis  of  moments  from 
the  roll  and  yaw  louver  inputs  approximately  rotates  with  collective  vector, 
an  analysis  was  made  whereby  the  rate  gyros  were  rotated  with  vector  also, 
such  that  gyros  and  controls  would  act  in  the  same  coordinate  system.  Fig¬ 
ures  77  and  78  show  the  effects  of  rate  gyro  rotation. 

In  Figure  77,  the  roll  gain  is  2,  and  it  is  seen  thatwN  seems 

N1 

not  at  all  affected  by  K  . 

y 


In  Figure  78  the  roll  gain  is  20,  and  both  u>  values  vary  over 
'  bout  the  same  range  for  the  yaw  gain  variation  from  2  to  40. 

An  open-loop  compensation  of  the  effects  described  above  can 
be  made,  similar  to  those  accomplished  in  the  mechanical  mixer  box. 

It  is  hoped  that,  if  the  above  effects  prove  objectionable,  the  problem  can 
be  minimized  by  simply  turning  the  yaw  SA  system  axis  off  when  0  be¬ 
comes  greater  than  15  to  20  degrees.  This  would  be  the  most  reliable 
modification  to  the  present  system. 

7.2  Analysis 

S  8 

Figures  79  and  80  are  plots  of  AC^  and  AC^  due  to  /?g 

and  0  inputs, 
v 

This  analysis  is  restricted  to  angle  of-attack  =  0,  which 

aligns  the  body  and  stability  axes.  Assume  trimmed  flight  prior  to  small 

incremental  louver  inputs,  /3  and  /3  . 

s  v 

Assuming  the  forces  due  to  louver  inputs  act  at  the  fan  hubs, 
the  resultant  rolling  and  yawing  moment  expressions,  in  terms  of  changes 
in  fan  lift  and  drag  are: 


Rolling  Moment  = 

<AVALrt>Yf 

(1) 

Yawing  Moment  =  Af  = 

(AV  AVYf 

(2) 

Where  =  wing  fan  moment  arm  from  fuselage  center  line,  L  =  lift, 
D  =  drag. 


It  is  necessary  to  determine  the  variations  of  lift  and  drag  on 
the  fans  due  to  louver  inputs. 

Using  Figures  79  and  80  we  define  A  L,  the  incremental 
change  in  lift,  as: 

AL  =  AC  SqSA  (3) 

Lt  i 

which  holds  for  both  wings. 


/ 
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For  example,  to  calculate  the  incremental  lift  change  on  the  right  wing: 

ALrt=*CLV7  <4> 

rt  2 


Going  one  step  further,  we  can  expand  the  expression  for  Ac  .  as: 

L  rt 


A  C. 


s 


B  C. 


rt  BP 


B  C 


s 


rt 


»'3s  srt 


(5) 
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The  same  expressions  hold  for  the  left  wing. 


For  a  positive  rolling  moment,  p  is  positive  and  p  is 

Srt  ^t 

negative,  and  they  are  assumed  equal  in  magnitude. 


For  a  positive  yawing  moment,  (3  increases  and  /? 

V£t  Vrt 


decreases. 


With  the  above  knowledge  at  hand,  it  is  possible  to  determine 
the  roll  and  yaw  characteristics  of  both  differential  stagger  and  differen¬ 
tial  vector  inputs  to  the  louvers. 


If  we  define  p  as  the  differential  stagger  variation  on  the 
right  wing,  and  p  as  the  differential  vector  variation  on  the  left  wing, 
a  positive  p  results  in  a  positive  rolling  moment,  and  a  positive  p 
results  in  a  positive  yawing  moment. 


For  a  positive  roll  input: 


AV*s>-  ALrt<<V 
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Knowing  P  =  -0  =  -  0 
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-2 


8  0 


(7) 


8C. 


s 


«<*.>-  VAf 


8/3 


8 


(8) 


8L(Pg, 

-  '  'V  Af 


8C, 


s 


8  0 
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8/3 
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(9) 


Now,  calculating  the  yawing  moment  effect  of  a  positive 


0  s  input 


'^8)  =  Yf(ADrt(/V  -AV  V 


(10) 


*</y  =  YA 

8  “2 


ac®  sc  8 

- 0  _ - 2—  p 

8Ps  ®rt  8  ^8  BJt 


(11) 


Knowing  0  =  -  0  =  “  0_ 

®LT  ®rt  8 


*(/*s)  =  YfqAf 


SC, 


s 


+  2  0 


s  80 


(12) 
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8C 


s 


•*(/y=  +Yfq  Af  P8 


D 


8/3 


(13) 


s 


•*('.>  _  aCD 

-  =  +  Y,q  A.  - 

8  /?  g  f  f  8/? 


s 


s 


=  >V 


(14) 


Similarly  for  positive  yaw  inputs: 


•*<V  =  Yf 


ADrt(^v>-AV^ 


(15) 


*<<>„)- VAf 


8C 


s 


D 


8C 


D 


dp  v  „ 
H  v  rt 


8^v  V 


(16) 


Knowing  fl  =  -  P  =“^v 
rt  it 


^(/Jv)  =  -YfqSAf^ 


8C 


8/3, 


(17) 


Sc  (  0^  =  Yf  A  Lit^  ^  v*  ’  A  Lrt^  V 


(18) 


(19) 


*</y  =  YA 

2 


8C, 


8  C 


s 


8  0, 


0 


vit  8*v  '  Vrt 


(20) 


*  (  0y)  =  Yfq  Af  /3 


8C, 


s 


8  P 


(21) 


(22) 


The  following  Table  is  a  listing  of 


8C, 


s 


0C 


Bp 


and 


D 


BP 


derived  from  Figures  79  and  80 ;  for  collective 


P  =  0,  10,  20,  30,  40  and  50° ,  and  P  trim  equal  to  27° 


s 


The  curve  of  Figure  79  was  used  to  determine  the  trim  T 


values  for  the  collective  p  values  given  above. 


P, 


s 


8, 


0 

10 

20 

30 

35 

50 


.98 

.955 

.90 

.855 

.8 


8ClM 


-.0093 


-.0120 


-.0115 


-.0086 


-.0050 


-.005 


8/ 


'd/  8Py 


BC 


s  j 


L/Wv 


BC 


.0048 


.0069 


.0048 


.0038 


.0038 


0 


-.003 


-.0118 


-.0130 


-.0138 


-.0140 


D/dpv 

-.0190 

-.0155 

-.0128 

-.0098 

-.0086 

-.0080 


(All  partials  per  degree) 


It  is  desired  to  determine  the  axis  of  moments  for  differential 
stagger  inputs  and  differential  vane  inputs. 

Figure  81  shows  the  axes  for  roll  and  yaw  moments.  A  posi¬ 
tive  rolling  moment  is  along  the  positive  "xM  axis,  and  a  positive  yawing 
moment  is  along  the  positive  "z"  axis. 
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The  resultant  of  a  rolling  and  a  yawing  moment,  •#,  has  a 


y2  2 

£  +  Jf  ,  and  is  oriented  in  the  xz  plane  at  an 

angle  6  -  tan-1  M/x  with  the  x  axis. 

Since  M  and  if  are  proportioned  to  the  — values,  the 

30 

angle  of  the  resultant  6  ,  may  be  computed  from  the  coefficients  alone. 


Thus 

if 


0. 


M 


0 


3C. 


30 

3CD 


s 

s 


30 


(23) 


(24) 


if 


0 


3C. 


(25) 


v 


ac 


D 


» 0 


(26) 


These  values,  as  well  as  the  resultant  magnitudes  and  angles  are: 


? 

V 

T  S 
c 

if 

* 

**  * 
V 

X 

Z  tan-1 

0 

1 

.0093 

0 

.0093 

z  o= 

0 

.019 

.019 

Z  90° 

10 

.98 

.0120 

.0048 

.0129 

Z  21.8  = 

-.008 

.0155 

.01745 

Z  117.3" 

20 

.90 

.0115 

.0069 

.0134 

Z  *1  = 

-.0118 

.0128 

.01735 

Z  132. 7” 

30 

•  9 

.  oosu 

.0048 

.00984 

Z  29.2  = 

-.0130 

.  0098 

.01625 

Z 142.9= 

35 

.  SO 

.  0050 

.  0038 

. 00627 

Z37.3  = 

0134* 

.  uOSfi 

.01025 

Z14^2; 

50 

.SO 

.  005 

.  00 3 H 

. 00627 

Z  37.3' 

-.01  ! 

.  OOS 

.0160 

Z  150 

To  convert  these  coefficients  to  rolling  and  yawing  moments, 
we  use  equations  9,  14,  18  and  22.  The  coefficients  must  be  multiplied 
by  Yf,sAf. 

g 

In  the  transition  region,  a  reasonable  value  for  q  is  300  psf, 
and  Yj  is  equal  to  5  feet.  is  42. 6  ft? 

Thus, 


YfqSAf  =  2  X  5  X  300  X  42. 6  =  63, 900  (27) 

For  the  case  of  collective  0  =  35° ,  the  value  used  in  the  coupling 

analysis,  we  can  calculate  the  values  of  £  ,  Jt  ,jC  ,  and  Jt  . 

P  s  *s  ^  v  ^v 

Thus,  at  0  =  35°,  0  =  27'  ; 

V  s 

£  =  . 005  X  63,  900  =  319.  5  ft  lb/deg  =  18, 300  ft  lb/rad  (28) 

^s 

Jt  p  =  .  0038  X  63,  900  =  243  ft  lb/deg  =  13, 925  ft  lb/rad  (29) 
s 

2  =  -.  0138  X  63, 900  =  -882  ft  lb/deg  =  -50, 500  ft  lb/rad  (30) 

P  v 

Jt  =  .  0086  X  63,  900  =  549. 5  ft  lb/deg  =  31,  486  ft  lb/rad  (31) 

H 

v 

Using  the  values  thus  derived,  we  can  now  develop  the  charac¬ 
teristics  for  an  inertial  two  degree- of-freedom  (roll  and  yaw)  system  with 
stability  augmentation. 

The  SA  system  will  be  connected  so  that  roll  rate  will  be  fed  into 
0  ,  and  yaw  rate  will  be  fed  into  0  A  180°  phase  shift  is  required  for 

stability. 

If  we  pick  a  compensation  network  with  the  transfer  function 

1  +  § — ,  and  let  K  =  roll  gain  and  K  =  yaw  gain,  then  we  can  write  the 
1  +  10S  R  y 

equations  of  motion  for  the  closed=loop  system.  Also,  since  the  charac¬ 
teristic  equation  of  the  system  is  the  desired  output,  we  can  write  the  force- 
free  equation  of  motion  (i.  e. ,  no  stick  inputs)  as  follows: 
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ps  = 


I 

x 


kr*/> 


-S. 


I 


rS  = 


I 

xz 


pS- 


S' 


/? 

s 


I 

z 


I 

z 


p 


+  r 


KY*e 

+ _  V 


I 

z 


0 


(34) 


(35) 


By  inspection,  and  knowing  that  the  characteristic  equation 
must  have  zero  coefficient  sign  changes  to  be  stable,  we  can  see  the  value 
of  the  constant  term  will  be 

[sVS  *  a  'Ws  % 

S  V  vs 


(1  +  sf  f 

(1  +  10S)* 


(36) 


For  stability,  the  value  of  this  term  must  be  greater  than  zero, 
so  the  primary  criterion  for  stability  is  that: 


Xa  X  n  * 


0_  0 


*  0 


s 


(37) 
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Substituting  in  the  constants  of  equation  28,  29,  30  and  31, 
it  is  seen  that  the  XV- 5 A  louver  system  at  least  satisfies  part  of  the 
requirements  for  a  stable  system. 


2 

Using  values  of  I  =  3930  slug  ft  ,  I  =  16, 180  slug 

2  X  Z 
I  =  830  slug  ft  ,  the  expansion  of  equations  34  and  35  result 
xz 

characteristic  equation: 


ft2,  and 
in  the 


98.  9S4  +  S3  (19. 78  +  96. 94  KR  +  25. 71  1^) 

+  S2  (28. 3Ky  +81.06  KyKR  +  106.  6kr  +  .989) 


+  S  ( 2.  57K  +  9. 69K  +  162K__K_)  +  81. 06K  K  =  0  (38) 

Y  K  Y  H  Y  xt 

Inserting  the  values  of  between  2  and  40  with  a  value  of 

K  =  2  in  equation  38  results  in  the  roots  plotted  in  Figure  75. 

R 

Figure  76  is  a  plot  of  the  same  roots  for  K  =  20. 

R 

As  seen  in  Figure  74,  and  as  before  noted,  the  angle  between 
the  /?  resultant  and  the  /?v  resultant  remains  approximately  90°  .  Thus 

analysis  of  the  effect  of  rotation  of  the  roll  and  rate  gyros  downward  with 

collective  / 3  was  made.  This  was  done  for  the  /?  =35°  case, 
v  v 

For  this  case,  the  rate  measured  by  the  roll  rate  gyro, 

is  equal  to:  p  =  p  cos  o  +  r  sin  o 
gyro  H  v  H  v 

Similarly,  the  rate  measured  by  the  yaw  rate  gyro  is: 

r  =  -p  sin  n  +  r  cos  r 
gyro  pv  H  v 

These  relationships  were  substituted  into  the  parts  of  equations 
32  and  33  involving  the  SA  system  feedback  terms.  The  characteristic 
equation  resulting  from  this  substitution  is: 

98.  92S4  +  S3  (19. 78  +  92.  2K_  +  164.  4K  ) 

R  Y 

(continued) 


Pgyro* 

(39) 

(40) 


/ 

'/ 
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+  s2  (.  989  +  101.  42Kr  +  .  03Kr2  +  80.  96^^  +  180. 8Ky) 


+  S  (9.  22Kr  +  .  06Kr  +  161. 3KyKR  +  16. 441^) 


+  ,03K_  +80.86  K  K  =  0 
R  Y  R 


(41) 


The  various  values  for  K  and  K  were  substituted  in  equation 

Y  R 

41,  and  the  results  are  plotted  in  Figures  77  and  78. 

For  a  drill,  a  35°  upward  gyro  rotation  was  investigated,  and 
this  resulted  in  an  unstable  system. 
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8.0  LIFT-FAN  PITCH  CONTROL  STRUCTURAL 
FEEDBACK  ANALYSIS 


8.1  Discussion 

When  stability  augmentation  and/or  autopilot  systems  are  used  in 
a  vehicle  which  cannot  be  assumed  inelastic,  control  signals  from  the 
sensors  can  cause  autopilot  loop  instability. 

The  stabilization  system  sensors  are  usually  oriented  to  detect 
rigid  body  motion,  and  their  output  phasing  is  based  upon  rigid  body 
considerations. 

Since  the  vehicle  body  is  a  flexible  structure,  the  sensors  in 
reality  measure  both  rigid  body  motions  and  elastic  body  deformations. 

If  the  elastic  deformations  have  sufficient  magnitude  and  the  wrong 
phasing,  control  loop  oscillation  can  result. 

The  means  used  to  describe  the  dynamics  of  the  flexible  body 
motion  are  similar  to  the  use  of  a  Fourier  Series  to  describe  an  arbi¬ 
trary  periodic  wave  form. 

The  motions  of  a  flexible  body  can  be  exactly  described  by  the 
sum  of  the  rigid  body  motions  plus  the  contribution  from  an  infinite 
number  of  "orthogonal  body  modes".  Each  of  these  "orthogonal"  or 
"normal"  modes  has  associated  with  it  a  "mode  shape"  and  "modal 
natural  frequency". 

One  property  of  these  "normal"  modes  is  that  the  kinetic  energy 
of  the  body  motion  will  be  the  sum  of  the  kinetic  energies  of  all  of  the 
normal  modes  plus  the  kinetic  energy  of  the  rigid  body  modes. 

If  the  flexible  body  is  excited  with  an  impulse  of  force,  the  re¬ 
sulting  body  vibrations  will  be  the  sum  of  contributions  from  the  various 
normal  modes. 

With  each  normal  mode  is  associated  a  "mode  shape",  and  also 
a  "generalized  mass"  for  the  mode.  The  mode  shape  is  a  plot  of  the 
modal  deflections  of  points  on  the  body,  normalized  to  any  convenient 


reference  deflection.  The  generalized  mass  for  a  mode  is  defined  as 
the  sum  of  the  differential  masses  of  the  body  times  the  square  of  the 
modal  deflections  at  the  various  mass  points  on  the  body. 


mt  =  2dm  [  <pW  ] 


(1) 


The  equation  of  motion  for  the  flexible  body  modes  are  given  as 
below,  in  Laplace  notation: 


F  <f>  <l)  +  M  a  (1) 
{  m 


(2) 


Where 


f  j  =  the  damping  of  the  ith  mdde 
u) .  =  the  natural  frequency  of  the  ith  mode 


q.  =  the  generalized  coordinate  of  the  ith  mode 


m.  =  the  generalized  mass  of  the  ith  mode 


i 

F  =  the  applied  force 


M  =  the  applied  moment 


<$> 


(i) 


(i) 


m 


=  the  normalized  deflection  of  the  ith  mode  at  the  point 
of  application  of  the  force,  F. 

=  the  slope  of  the  normalized  deflection  curve  of  the  ith 
mode  at  the  point  of  application  of  the  moment,  M. 


In  addition,  2 
force". 


F* ,(i>  ♦  Mo  « 
f  m 


is  defined  as  the  "generalized 


In  general,  in  autopilot  analysis  we  are  interested  in  solving  for 
the  generalized  coordinates,  because  the  sum  of  the  generalized  coor¬ 
dinates  of  a  point  on  the  body  will  determine  the  deflection  of  that  point 
on  the  body. 

We  can  show  in  block  diagram  form  the  effect  of  the  rigid  + 
flexible  body  modes. 


Let  us  assume  that  we  are  interested  in  the  effect  of  the  body 
flexibility  upon  the  pitch  rate  gyro.  Further  assume  that  the  pitch  rate 
gyro  signal  is  finally  fed  to  the  nose-fan  thrust  reverser.  This  is  shown 
in  Figure  82. 

Thus,  the  pitch  rate  gyro,  located  at  station  Xg,  measures  the 
total  pitch  rate  at  this  station.  This  signal  is  then  fed  to  the  nose-fan 
thrust  reverser  servo  located  at  approx.  X^,  which  modifies  the  body 
motion  in  the  desired  (it  is  hoped)  manner. 

This  is  shown  in  block  diagram  form  in  Figure  83. 

9g  is  defined  as  the  total  pitch  angle  at  the  gyro  station,  and  & 
is  the  nose-fan  thrust  reverser  angle,  measured  from  full  closed.  ^ 

In  general,  an  infinite  number  of  flexible  body  modes  are  not 
required  to  represent  the  elastic  body. 

In  our  case,  the  first  8  modes  were  investigated.  The  effect  of 
the  7th  &  8th  modes  was  found  to  be  negligible. 

Referring  to  Figure  83,  it  is  seen  how  the  total  pitch  angle  at  the 
gyro  station  is  computed. 

In  the  case  of  the  XV-5A,  the  excitation  for  the  symmetrical  flex¬ 
ible  body  modes  is  provided  by  the  nose-fan  force.  This  force  consists 
of  two  components.  The  first  is  a  force  due  to  thrust  modulation  of  the 
nose-fan  by  means  of  thrust  reverser  door  movements.  The  second  is 
a  reaction  force  due  to  the  acceleration  of  the  thrust  reverser  doors . 
Figure  84  shows  the  rigging  of  the  nose-fan  doors. 

The  door  actuating  mechanism  is  such  that  for  simusoidal  inputs 
to  the  doors,  the  reaction  due  to  door  acceleration  becomes  equal  to  the 
change  in  thrust  due  to  door  displacement  at  a  frequency  of  about  5  cps. 

Above  this  frequency  the  reaction  due  to  door  acceleration  be¬ 
comes  dominant. 

Due  to  the  iayout  of  the  door  actuation  mechanism,  essentially 
zero  moment  is  applied  at  the  nose-fan  station  due  to  door  operation. 

The  generalized  force  used  to  determine  the  flexible  body  coordinates, 
q[t  is  thus  2  F  <p  NJ1  \ 

P 

Once  the  q^'s  have  been  determined,  their  contributions  to  the 

slope  (pitch  angle)  at  the  gyro  station  are  found  by  multiplying  the 
values  by  the  slope  of  the  <f>^  curves  at  the  gyro  station,  . 


The  total  pitch  angle  at  the  gyro  station  is  then  the  sum  of  the 
rigid  body  pitch  angle  plus  the  contributions  from  the  flexible  body  modes. 

The  values  are  functions  of  S,  as  is  shown  in  equation  (3) 


F  <p 


(i) 


Mo 

+  m 


(i) 


m^  (S2  +  2ttw  S4-W  ^ 


(3) 


We  can  then  write  the  total  pitch  angle,  £  9  ,  at  the  gyro  sta¬ 

tion  due  to  the  flexible  body  modes  and  add  them  to  the  rigid  body  0»  as 
is  shown  in  equation  (4) . 


9  =  9  ,  _  +  2 

g  rigid 


(»'**"'.*)*  a* 

n^  (S^  +  2^ituiS4  ^ 


or,  in  our  case,  assuming  the  rigid  body  is  inertial,  we  write: 


(4) 


9 

g 


M. 


+  £ 


5  <h 

P-NF 


I  S 

y 


2  2 
m^(S  +2£  ) 


Where: 


(X-X„  _ 
v  NF  C.G. 


) 


(5) 


(6) 


and  F  is  the  force  at  the  nose  fan  station  due  to  <5  ,  and  Includes  the 

*p  p 

effect  of  nose  fan  door  inertia. 


As  is  seen,  equation  (S)  is  the  sum  of  a  series  of  fractions,  and 
to  get  the  9  /  6  transfer  function  we  must  determine  the  sum  of  these 

8  p 

fractions.  The  lowest  common  denominator  of  these  fractions  is  the 
product  of  all  the  fraction  denominators,  since  they  are  all  different. 
For  each  gyro  location  investigated,  X_  ,  the  numerator  terms  vary, 

but  the  denominator  terms  remain  the  same. 


When  this  operation  is  performed,  a  numerator  term  results,  and 
this  numerator  term  is  a  function  of  S,  the  Laplace  operator.  The  numer¬ 
ator  is  then  generally  a  polynomial  in  S,  and  as  such  has  roots  which 
result  in  the  polynomial  having  zero  value.  These  roots  are  the  zeros 
of  the  8g/d  transfer  function,  while  the  roots  of  the  denominator  cause 
the  transfef*  function  to  have  an  infinite  value,  and  are  called  the  poles 
of  the  transfer  function.  As  previously  stated,  the  zeros  vary  with  gyro 
location,  but  the  poles  are  invariant. 

In  general,  the  0  /b  transfer  function  is  favorable  if  the  zeros 
of  the  transfer  function  <?cci?r  at  lower  frequencies  than  their  associated 
poles.  This  is  to  say,  the  roots  of  the  transfer  function  should  be  aligned 
zero,  pole,  zero,  pole,  etc. ,  as  S  increases.  This  configuration  of  poles 
and  zeros  results  in  a  phase  characteristic  which  tends  to  lead  the  rigid 
body  phase  angle. 

The  local  pitch  angle/nose  fan  force  response  has  been  calculated 
for  gyros  located  at  Sta.  90,  Sta.  120,  and  Sta.  150. 

Bode  plots  for  these  responses  are  shown  in  Figures  85,  86  and  87. 

Figure  85  is  the  plot  for  the  Sta.  90  gyro  location.  This  configura¬ 
tion  results  in  a  large  additional  phase  lag  in  the  100-150  rad/sec  region. 

Figure  86  is  for  the  Station  120  gyro  location.  It  is  seen  that  the 
phase  angle  never  becomes  greater  than  180°  lagging,  which  is  the  phase 
angle  of  the  assumed  rigid  body  mode. 

Figure  87  is  the  plot  for  the  Station  150  gyro  location.  This  loca¬ 
tion  is  in  the  equipment  compartment  area,  and  was  the  originally  selected 
gyro  location.  This  configuration  also  has  large  phase  lags  around  100 
rad/sec. 

In  order  to  combine  the  flexible  body  with  the  SA  system,  several 
of  the  parameters  were  assigned  reasonable  values. 

The  rate  gyro  was  given  a  natural  frequency  of  26  cps,  with  a  damp¬ 
ing  ratio  of  0.3  critical. 

The  servo  was  given  a  third-order  response  with  the  first  order 
time  constant  of  1/20  second,  and  the  second  order  having  a  natural  fre¬ 
quency  of  20  cps  and  damping  ratio  of  0.3. 
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The  nose  wags  body  transfer  function  is  a  double  zero  at  5  cps, 
with  a  gain  constant  of  -1450. 

The  gyro  compensation  network  has  the  transfer  function 

Amplifier  input  1  +  S 

Gyro  output  =  1  +  10S 


This  corresponds  to  a  ratio  of  .  1,  which  is  the  worst  case  as  far 
as  structural  feedback  is  concerned. 


Figure  88  is  a  Bode  plot  of  the  response  of  the  servo  and  nose  fan 

thrust  reverser  response,  (F/F  ) , 

o  Ik  _ 


Figure  89  is  a  Bode  plot  of  the  complete  system  for  the  Station 


120  gyro  location.  In  this  case  the  allowable  gain,  $ 
deg/sec. 


p/0 


is  13  degrees/ 


g 


Figure  90  is  the  same  as  Figure  89,  with  one  exception.  In  this 
case,  a  parallel  MT"  notch  network  has  been  added  with  a  null  frequency 
of  100  radians/sec.  This  network  has  the  transfer  function: 


1 


2S 


+ 


S 


2 


co 


2 

o 


In  this  case,  the  allowable  gain  can  be  raised  to  70,  or  37  DB. 

Figure  91  is  a  Bode  plot  of  the  complete  system  less  the  notch 
network  for  the  Sta.  150  gyro  location.  The  maximum  gain  attainable 
with  this  configuration  is  2. 9. 


Figure  92  is  a  Bode  plot  for  the  Sta.  150  configuration  with  the 
addition  of  the  100  rad/sec.  notch  network.  Here  a  maximum  gain  of 
40,  or  32DB  is  attainable. 


Figures  93,  94,  95  and  96  show  the  modal  shapes  along  the  fuse¬ 
lage  centerline  for  the  first  8  body  modes.  It  is  seen  that  Station  120 
is  an  advantageous  gyro  location  because  the  slopes  of  the  first  two  modes 
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are  small,  the  slopes  of  the  third  and  fourth  modes  are  small,  the  slopes 
of  the  fifth  and  sixth  modes  are  zero,  and  the  slopes  of  the  seventh  and 
eighth  modes  are  negligible. 

For  this  work  the  rigid  body  was  assumed  to  be  purely  inertial. 
This  is  not  true  at  low  frequencies,  but  at  the  frequencies  of  interest  in 
this  investigation  (10  cps  and  up)  this  is  a  valid  assumption. 

In  general,  the  vehicle  will  be  flown  at  various  gain  and  ratio 
settings,  and  the  objective  of  this  analysis  is  to  ensure  that  within  the 
gain  capabilities  of  the  SA  System  the  rigid  body  modes  will  be  the  limit¬ 
ing  factor  on  system  gain.  A  system  gain  of  about  100  is  attainable  from 
the  stabilization  package,  but  the  short  period  will  go  unstable  at  a  gain 
of  around  50.  It  is  thus  desirable  to  be  able  to  attain  a  gain  of  over  50 
without  structural  feedback  problems. 

In  addition  to  relocating  the  gyros  to  Station  120,  the  incorpora¬ 
tion  of  the  parallel  "T"  notch  network  has  been  recommended.  This  net¬ 
work  is  light  and  simple,  and  affords  zero  transmission  at  its  design 
frequency.  It  is  useful  to  reduce  the  system  amplitude  at  such  a  point 
where  a  troublesome  resonant  peak  occurs. 

Preliminary  tests  on  the  unloaded  nose  fan  door  servo  indicate 
that  the  damping  of  the  second-order  portion  of  the  response  might  be 
higher  than  previously  estimated.  The  addition  of  the  servo  inertia  load 
should  both  lower  the  frequency  of  the  second-order  and  the  damping. 

Figure  97  is  a  plot  of  phase  and  amplitude  for:  the  servo, 

w  =120,  £  =  .3,  T  =  1/20  sec. ;  the  rate  gyro,  w  =  163,  £  =.3;  the 
n  n 

compensation  network,  t  numerator  =  1  sec. ,  r  denom  10  sec. ;  and 
the  N.  W.  B.  zero,  w  N  =  31. 5,  £  =  0. 0. 

Figure  98  is  a  phase-amplitude  plot  for  the  same  configuration 
as  Figure  97,  plus  the  addition  of  a  parallel  "T"  notch  network  with  the 
null  at  100  rad. /sec. 

Figures  99,  100  and  101  give  the  complete  system  block  diagram 

for  the  Sta.  90,  120  and  150  gyro  locations.  Due  to  the  large  number  of 

second-order  terms  in  the  6  / 6  transfer  function,  a  shorthand  notation 

g  p  . 

is  used.  The  expression  9o  =  1/  w  =A  (3)  is  defined  as  equivalent 

9.  /  fc=B 
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to  the  expression 


8. 2  Analysis 

The  inertia  of  the  nose  fan  doors  has  considerable  bearing  on  the 
structural  feedback  problem.  This  inertia  is  reflected  back  through  two 
bellcranks  to  the  nose  fan  door  actuator,  which  reacts  vertically  against 
the  aircraft  structure. 

The  inertia  of  a  nose  fan  door  must  first  be  calculated  referred 
to  the  door  hinge  axis. 

According  to  the  latest  weights  data  available,  the  freight  of  one 
door  is  22  lbs. ,  I  of  a  single  door  is  equal  to  .  115  slug  ft  ,  the  c.  g.  of 
the  door  is  at  X  =*56. 8,  ~Z  =  81. 9,  and  Y*  =  16. 2.  The  front  hinge  is 
located  at!T=  34.  55,  T=  12. 0,  and“z  =  83. 16. 

The  aft  hinge  is  located  at  X  =  79. 8,  Y  =  17.34,  "z  =  77. 82. 

To  transfer  this  inertia  to  the  door  hinge  line,  we  must  deter¬ 
mine  the  distance  from  the  door  c.g.  to  the  hinge  line.  This  is  done 
most  easily  by  means  of  descriptive  geometry,  as  in  Figure  102. 

Frorr  Figure  102,  the  distance  between  the  door  c.g.  and  door 
hinge  line  is  l  inches.  We  can  thus  calculate  the  door  moment  of  inertia 
about  the  door  hinge  line.* 

9 

(9  V  /  22  2  \  2 

W!  (32^2)  =  •115  +  *0191  -  -134  slue ft  (i) 

We  must  now  calculate  the  effective  inertia  reflected  to  the  nose 
fan  door  actuator,  and  refer  it  to  6  ,  the  door  position  angle. 

p 

Figure  64  is  a  1/4  scale  layout  of  the  door  rigging. 

A  door  position  of  33°  from  full  thrust  was  picked  as  nominal  for 
this  investigation,  and  this  is  the  position  used  in  Figure  84. 

As  shown  on  Figure  84,  the  force  required  from  the  actuator  for 

o  V 

nose  fan  door  acceleration  at  5p  =  33  is  -5.36  o  p  I  door*  *or  eac^ 


door.  The  total  is  equal  to  -10. 72  6  ID0OR-  This  force  is  vertical* 
and  an  upward  actuator  force  results  in  a  positive  6  p.  The  reaction  on 
the  airplane  is  down,  which  is  a  negative  force  for  a  positive  6  p. 

In  addition,  the  reactions  on  the  two  hinge  points  amounts  to 

-1. 34  5  pI^oOR  f°r  eac^  slc*e  the  airplane,  for  -2. 68  6  p^DOOR  total . 

•  •  •  • 

Thus,  the  total  reaction  force  due  to  6  is  equal  to  -13.46  pi^QQ^  or 
“kl  «p 

These  figures  assumed  the  door  c.g.  was  on  the  hinge  line, 
whereas  it  is  actually  2"  displaced.  This  effect  is  negligible  on  the 
results. 

We  have  linearized  the  thrust  change  due  to  6p  such  that  75°  6  p 
results  in  -1700  lbs.  thrust  change,  or  -1300  lbs/rad. 

We  can  thus  express  the  total  force  from  nose 

Fibs.  =  1300  6  -  1.86 

P  P 

•  • 

or  Fibs.  =  -1300  (6  +7=sjr) 

'  p  722  ' 

expressing  (3)  in  Laplace  notation, 

2 

-1300  (l+|^)  lbs/rad. 

F  S2 

-  =  -1300  (1  + - -)  lbs/rad. 

6p  26. 9 

The  structural  feedback  analysis  used  earlier  data,  and  as  a  re¬ 
sult  the  transfer  function  used  was: 

s2 

F/6  =  -1450  (1  + - -)  lbs/rad.  (6) 

P  31.5 


(4) 

(5) 


fan  door  deflection: 

(2) 

(3) 
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Flgurt  83  Block  Diagram  o / 
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Figure  84  Nose  Fan  Door  Rigging  -  One  Side 
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FUSELAGE  STATION,  INCHES 

Figure  93  Normalized  Deflection  vs  Fuselage  Station  for  1st  and  2nd  Free-Free  Body  Mode 
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Figure  94  Normalized  Deflection  vs  Fuselage  Station  for  3rd  and  4th  Free-Free  Body  Mode 


STH  MODE,  f-  16.6  Op*,  »  •  1.44 
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MODE,  f  •  30.4  op*.  ■  1.75 


FUSELAGE  STATION.  INCHES 

Figure  95  Normalized  Deflection  vs  Fuselage  Station  /or  5t/>  and  6th  Free-Free  Body  Mo  dot 
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/ 
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7TH  MODE,  f-  23.9  cps,  m?  -  .5*08 
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Figaro  96  Normalized  Deflection  vs  Fusolago  Station  /or  7th  and  8th  Free-Free  Body  Modes 
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9.0  YAW-ROLL- PITCH  EULER  ANGULAR  AND  DEFLOREZ 

DISPLAY  RATES  DERIVATION 

9.1  Discussion 


Early  in  the  XV-  5A  program,  the  need  became  apparent  for 
a  more  realistic  visual  pilot  display  if  the  simulation  program  was  to 
have  meaningful  results. 

A  nationwide  survey  trip  was  made  to  investigate  various 
display  concepts,  including  both  fixed  and  moving  base  presentations. 

The  most  promising  system  was  the  point  light  source  system  developed 
by  the  DeFlorez  Company  of  Englewood  Cliffs,  New  Jersey,  a  fixed 
base  system. 

In  this  system,  a  transparency  containing  the  various  terrain 
features  is  moved  about  under  the  light  source,  and  the  point  source 
projects  through  the  transparency,  resulting  in  a  realistic  display  on 
the  spherical  screen  surrounding  the  simulator  cockpit.  If  the  pilot's 
eyes  were  at  the  location  of  the  point  source,  an  exact  representation 
would  be  afforded.  Depending  on  the  radius  of  the  screen,  some  distance 
is  allowable  between  the  eyes  of  the  pilot  and  the  point  light  source. 
Translation  of  the  transparency  simulates  vehicle  translation,  and  ro¬ 
tation  of  the  transparency  simulates  vehicle  rotations. 

A  specification  for  a  system  such  as  this  containing  the 
various  dynamic  response  requirements  was  prepared  by  XV- 5A  per¬ 
sonnel,  and  the  device  was  constructed  by  the  DeFlorez  Company. 

The  transparency  is  six  feet  square,  and  the  active  region  in 
X  and  Y  is  *  2  feet.  The  Z  travel  of  the  transparency  with  respect  to 
the  light  source  is  5  inches  total.  The  active  region  in  rol1.  is  approx¬ 
imately  -  25° ,  and  *17°  in  pitch.  Yaw  travel  is  unlimited. 

It  is  thus  possible  to  "scale"  the  linear  motions  of  the  device 
to  any  desired  ratio,  depending  upon  the  vehicle  requirements. 

One  of  the  limitations  of  the  scale  factor  chosen  is  the  max¬ 
imum  velocity  of  the  vehicle.  If  it  is  desired  to  be  able  to  fly  continuously 
around  the  transparency,  the  radius  of  turn  for  the  maximum  allowable 
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display  devloe  bank  angle  must  be  considered, 
the  radius  of  turn  Is  given  In  eq.  1. 

V2 

R  "  g  tan  $ 


For  a  coordinated  turn, 


(1) 


where  R  is  radius  of  turn  in  feet,  V  Is  velocity  in  ft/sec. ,  g  is  32. 2 
ft/sec2,  and  <p  Is  the  roll  angle. 


For  the  display  uevice,  <p  is  25° ,  so  equation  (1)  becomes  equa¬ 


tion  (2). 


(2) 


Now,  we  know  that  the  maximum  circle  we  can  put  in  a  4  foot  square 
has  a  2  foot  radius. 


If  we  pick  a  scale  factor  K, 
is  the  scaled  equation: 


then  we  can  write  equation  (3)  which 


KR  = 


(3) 


or 


or 


V2  =  K.  15*  2’ 
max 


V  2 
max 

30 


30K. 


(4) 

(5) 


If  we  thus  know  the  maximum  velocity  in  a  given  flight  regime,  we 
can  determine  minimum  scale  factor  which  allows  continuous  flight. 


In  general,  it  would  not  be  desirable  to  be  required  to  fly  contin¬ 
uously  in  a  circle,  so  a  good  value  for  would  allow  a  circle  to 

be  made  in  1/2  of  the  transparency  action  range,  or  1  foot  in  radius. 
Then  becomes,  in  equation  (6): 


V  2 
max 


(6) 


15 


There  is  another  important  requirement  limiting  the  settle 
factor,  and  this  is  the  required  minimum  altitude  of  flight,  hi  a  VTOL 
vehicle,  in  the  low- speed  and  landing  regime,  it  Is  desirable  to  be  able 
to  start  from  ground  level,  and  fly  the  complete  low- speed  envelope. 

In  the  case  of  the  XV-5A,  the  low-speed  envelope  includes 
speeds  up  to  conversion  speed,  or  around  100  knots.  We  may  thus 
determine  K  win  for  100  Kts  TAS.  If  we  are  to  be  able  to  land,  we  must 
be  able  to  attain  a  minimum  altitude  equal  to  the  level  of  the  pilot's  eye 
when  the  airplane  is  on  the  ground.  In  the  XV- 5A  the  pilot's  eye  level 
is  approximately  7  feet. 


The  limit  on  display  device  minimum  altitude  is  a  result  of 
the  fact  that  the  point  light  source  has  a  finite  diameter,  which  at  close 
proximity  to  the  transparency  causes  a  loss  of  resolution  of  the  project¬ 
ed  image  (fuzziness).  This  limit  has  been  set  at .  08  inches.  We  may 
thus  determine  Kmax  such  that  this  value  corresponds  to  7' or  84  inches. 
Thus,  in  equation  (7): 


.08" 


(  K 

\  max 


84" 


K 

max 


8400 

8 


1050:1 


(7) 


Knowing  V  =  170  fps,  we  calculate 
max 

K  =  =  1900:1 

min  i5 

In  the  hovering  and  transition  regime,  then,  we  have  two 
conflicting  requirements  on  scale  factor. 


We  thus  compromised  by  choosing  a  scale  factor  of  1500:1, 
which  allows  a  10'  minimum  altitude  and  still  allows  continuous  flight, 
although  somewhat  cramped. 

A  scale  factor  of  250:1  was  chosen  for  the  hovering  regime, 
and  50, 000:1  for  conventional  flight. 

Angles  are  reproduced  exactly,  so  no  conversion  is  necessary. 
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Vehicle  motions  are  thus  represented  by  motions  of  the  pro¬ 
jection  of  the  point  light  source  through  the  transparency  and  onto  the 
screen. 


The  earth  frame  of  reference  is  the  transparency,  and  the 
frame  of  the  device  (along  with  the  simulator  cockpit)  becomes  the  ve¬ 
hicle.  We  then  have  motions  of  the  earth  coordinates  which  look  like 
vehicle  motions. 

The  projector  must  be  driven  by  the  proper  computer  inputs 
to  truly  represent  the  vehicle  motions.  The  angles  through  which  the 
projector  must  rotate  are  referenced  to  the  earth  coordinate  system, 
and  are  thus  Euler  angles. 

In  the  projector,  the  pitch  axis  is  rigidly  attached  to  the 
support  structure.  The  roll  axis  is  carried  by  the  pitch  system,  and 
yaw  is  attached  to  roll. 

If  we  start  at  the  transparency,  which  actually  represents 
the  earth,  the  first  rotation  is  yaw,  the  second  is  roll,  and  the  last  is 
pitch. 


Thus,  we  have  a  set  of  Yaw,  Roll,  Pitch  Euler  angles,  and 
it  is  necessary  to  compute  the  rate-of-  change  of  these  Euler  angles  in 
terms  of  the  Euler  angles  and  the  body  angular  rates,  since  body  angular 
rates  are  computed  on  the  analog  computer. 

In  addition,  the  Euler  angles  must  be  used  to  calculate  ground 
velocities  in  terms  of  body  velocities. 

Section  9. 2  following  shows  a  derivation  of  the  Yaw-Roll- Pitch 
Euler  angles  and  direction  cosines,  and  the  derivation  of  the  ground  ve¬ 
locities  from  body  velocities. 

These  equations  are  used  to  generate  the  Euler  angles  required 

to  drive  the  projection  device;  as  well  as  to  geLerate  the  S  ,  S  and  S 

x  y  z 

components  for  drive  of  the  linear  degrees  of  freedom. 


9.2 


Analysis 


The  mechanical  layout  of  the  De  Florez  projection  device 
requires  the  use  of  Yaw,  Roll,  Pitch  direction  cosines  and  Euler  angles. 

The  derivation  of  the  direction  cosine  matrix,  Euler  angular 
rates,  and  inertial  velocities  for  this  system  of  Euler  angles  is  given 
below. 


We  use  a  right  hand  coordinate  system,  where  rotations  in  the 
order  of  X,  Y,  Z  are  considered  positive.  That  is,  a  rotation  of  X  into 
Y,  Y  into  Z,  or  Z  into  X  is  a  positive  rotation. 

For  the  Yaw  rotation,  referring  to  Figure  103,  104  and  105, 
we  can  write, 


X„  =  X  cos  ^  +  Y  Sin  * 
B1  e  e 


Y  =-X  sin  f  +  Y  Cos  * 
6  6 


z„  =  z 

B1  e 


For  the  Roll  rotation 


X„  =  X„ 

B  B 
2  1 


Y  =  Y  Cos  tp  +  Z  Sin  <p 
2  1  1 


Z  =  -Y  Sin  <#>  +  Z  Cos  <f> 

B2  B1  B1 


and  for  the  Pitch  rotation, 


X  =  X  Cos  0  -  Z  Sin  0 

B_  Brt  ii 

3  2  2 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 


0 
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(8) 


Y_  *  Y 
B  B 
3  2 


Z  ■  X _Sln  9  +  Z  Cob  # 
B3  B2  B2 


(9) 


We  can  write  these  equations  in  matrix  form  as  below: 


Cos  t 

Sin  ^  0 

X 

e 

= 

-Sin 

Cos  f  0 

Y 

e 

(10) 

1 - 

tSJ 

0 

M 

1 _ 

0 

0  1 

Z 

e 

V 

1 

0  0 

Yb2 

= 

0 

Cos  4  Sin  0 

\ 

(11) 

Zb2 

0 

-Sin  0  Cos^ 

- 1 

X 

0 

CO 

J 

Cos  # 

0  -  Sin# 

V 

YB3 

= 

0 

1  0 

(12) 

% 

Sin  # 

0  Cos# 

0 

to 

1 _ 

Combining  the  above,  and  solving  for  X  ,  Y  and  Z  in  terms  of 

U3  3  a 

X  ,  Y  ,  Z  ,  there  obtains: 
e  e  e 
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(13) 


I 

I 

i 

K 

f 

i 

i 

i 

i 


f 


B, 


B. 


B. 


Cos  e  0  -Sin  e 


1  o 


Sin  6  0  Cose 


10  0 


0  Cos0  Sin<f> 


0  -Sin0Cos0 


or 


B. 


B„ 


B. 


<p 


0 


X 


Expanding  the  product  6  <p  0,  we  get: 


Cos  0  Sin  0 


-Sin  0  Cos  00 


0  0  1 


X 


(14) 


V 

^ Cose  Cosip -  SinfiSin0  Sin  0^  Cos  $  Sin  0  +  Sin0  Sin0  Cos^j  (-SinflCos#) 

X 

e 

\ 

- 

|-Cos0Sin0^  (co80  Cos 0^  (sin  0  ) 

Y 

e 

V, 

(sin  0  Cos  0  +Sin  0 Cos 6  Sin 0 )  (sin 9  Sin0  -Cos  6  Cos  0 Sin  0  j  (cos6Cos0^ 

Z 

0 

m 

«.  _ 

(iS> 


I 

I 

I 


Or,  with  the  following  definitions  - 


The  various  t ,  m,  and  n  values  are  the  direction  cosines  of  the 
angles  between  the  body  axes  and  the  inertial  axes. 

Therefore,  for  example,  l  ,  m  ,  and  n  are  the  cosines  of  the 

X  X  X 

angles  between  the  body  X,  Y,  and  Z  axes  and  the  inertial  X  axis. 

The  direction  cosine  matrix  is  orthogonal,  and  thus  the  inverse  is 
equal  to  the  transposed  matrix,  and  the  determinant  of  the  matrix  has  a 
value  of  ±  l. 

We  can  then  modify  equation  (16)  to  solve  for  the  earth  components 
in  terms  of  the  body  components. 


X„  ,  Y_  and  Z„  are  the  body  coordinates  after  the  three  revolutions. 
B3  ®3  B3 

We  can  also  use  the  direction  cosine  matrix  to  determine  ground  rates 

from  body  velocities.  If  8  ,  S  and  S  are  the  ground  velocities  along  the 

x  y  z 

inertial  X,  Y  and  Z  axes,  and  u,  v,  and  w  are  the  body  velocities  along  the 
body  X,  Y  and  Z  axis,  then: 
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(18) 


m 

x 


m 

y 


m 


z 


■■ 

M  — 

n 

u 

X 

n 

V 

y 

n 

w 

z 

_  _ 

In  order  to  compute  the  Euler  angles  for  use  with  the  display,  it  is 
necessary  to  derive  the  Euler  angular  rates  and  then  Integrate. 


If  we  denote  the  matrix  of  direction  cosines  In  equlv.  (18)  as  X  we 
can  write  an  expression  for  the  rate  of  change  of  this  matrix  as  follows: 


(19) 


This  is  so  because  we  know  the  determinant  of  the  matrix  has  the 
value  -  1.  Thus  the  only  change  possible  in  the  matrix  must  be  due  to 
a  rotation. 


If  we  let  the  body  unit  vectors  equal  X,  Y  and  Z,  we  thus  can  write: 

A+wX  A  =  0  (20) 

or  A  =  -  u)X  A  (21) 

and  then,  ifw  =  p,  co  =  q,  and  =  r,  (22) 

A_  Y_ 

B  B  B 


~  A  * 

A  =  X  i 


A  .  A. 

+  Ym  +  Zn 
z  z  z 

;  a 

Xi  ■  rm  -qn 
z  z  z 

m  =pn  -r  Jo 
z  z  z 


A 

■  « 

A 

■  ■ 

A  ‘ 

=  -X 

qn  -rm 
z  z 

•  m 

-Y 

r  t  -pn 
z  z 

-z 

pm  -q  £ 
z  z 


(23) 

(24) 


•  * 

Vq  ^  z  pmz  (and  ditto  for  the  X  and  Y  values)  (25) 


using  the  values  for  £  ,  m  and  n  from  equations  (15)  and 

z  z  z 

(16)  we  write: 

:  d 

i  =rm  -qn  =  —  (-Sin0Cos0)  =  rSin0-qCos0Cos  0  (26) 

z  z  z  at 

m  =pn  -r  Z  -  -$j-  (Sin0)  =  pCos  0  +rSin  0  Cos  0  (27) 

z  2  z  at 

•  d  A 

n  =q£  -pm  =  -r*  (Cos  0  Cos0)  =  -q  Sin  0  Cos0-  pSin0  (28) 

Z  Z  Z  Gt 

Performing  the  derivative  operation 


rSin0  -q  Cos0Cos0  =  -  0  Cos0  Cos  0+  <f>  Sin0  Sin</>  (29) 
pCos  8  Cos  0  +  r  Sin  8  Cos  0  =  0  Cos  0  (30) 

-qJSin  8  Cos  0  -p  Sin  0  =  -  8  Sin  8  Cos  <p-  <f>  Cos  8  Sin  0  (31) 


from  equation  30,  we  see: 

0  =  p  Cos  8  +  r  Sin  8 
Using  equation  (29),  and  solving  for  8  : 


(32) 


^  q  Cos  8  Cos0-rSin0  +0  Sin0  Sin0  /oov 

Cos  0  Cos  0  (33) 


*  r  tan  <6 

0  =  q  -  — - 7—+  0  tan  0  tan  0 

Cos  0 


Substituting  0  =  p  Cos  0  +  r  Sin  0  ; 


(34) 


0  =  q  -  r  tan  0  Sec  0  +  p  tan  0  Sin  0  +  r  Sin  0  tan  0  tan  0  (35) 

Sin2  0 -1 


$  =  q  +  p  tan  0  Sin  0  +  r  tan  0 


Cos  0 


0  =  q  +  (p  Sin  0  -  r  Cos  0 )  tan  0 


Now,  we  can  write 


(36) 

(37) 


mw  =  p  n  -r-0 

y  y  y 


(38) 
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m  =  Cob  0  Cos  ip 

y 


(39) 


-37  (Cos  0  Cos  Ip  )  = 
at 

p  (Sin  9  Sin  ip  -  Cos  9  Cos  ip  Sin  <p ) 

-r(Cos  6  Sin  ip  +  Sin  9  Sin  <p  Cos  ip )  (40) 

•  • 

-  <p  Sin  <p  Cos  4>  -  ip  Cos  <p  Sin  $  = 

p  (Sin  9  Sin  4>  -Cos  9  Cos  ip  Sin  <p ) 

-r(Cos  9  Sin  ip  +  Sin  6  Sin  <p  Cos  ip  )  (41) 

t 

and  this  is  solved  for  ip  to  give: 

^  =  Sec  0  (r  Cos  9  -p  Sin  9  )  (42) 

The  Euler  angular  rates  are  then  given  in  equations  (32,  (37), 


and  (42),  repeated  below 

0  =  p  Cos  9  +  r  Sin  9  (43) 

9  =  q  +  (p  Sin  9  -r  Cos  9  )  tan  <p  (44) 

ip  =  Sec  <p  (r  Cos  9  -p  Sin  6  )  (45) 

and  we  can  also  write  the  direction  cosines: 

£  -  Cos  6  Cos  ip  -  Sin  9  Sin<J f>  Sin*P  (46) 

X 

m  =  -Cos  <p  Sin  *  (47) 

X 

n  =  Sin  9  Cos  ip  +  Sin  ip  Cos  9  Sin  0  (48) 

X 

£  =  Cos  9  Sin  ip  +  Sin  8  Sin  0  Cos  0  (49) 

y 
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m 

y 

=  Cos  <p  Cos  ip 

(60) 

n 

y 

-  Sind  Sin  ip  -  Coed  Cos  ip  Sin  <p 

(61) 

i 

z 

=  -Sin  0  Cos  <p 

(52) 

m 

z 

—  Sin  <p 

(53) 

n 

* 

=  Cos  0  Cos  <p 

(54) 

Due  to  the  limited  rotation  of  the  visual  display  in  pitch  and 
roll,  it  is  not  necessary  to  compute  the  direction  cosines  exactly. 

The  following  approximations  are  used  when  driving  the 
visual  display: 

$  =  p  +  r  8  (55) 

0  =  q  +  p  0  <f>  -  r  <p  (56) 

t  j*„n  Q 

tp  =  where  all  angles  are  in  radians  (57) 

and 


=  Cos  ip  -  8  <p  Sin  ip 

(58) 

=  Cos  <p  Sin  ip 

(59) 

«0Cos  ip  +  0  Sin  ip 

(60) 

■  Sin  ip  +  6  <p  Cos  ip 

(61) 

m  COS  Ip  COS  <p 

(62) 

«  8  Sin  ip  m<p  Cos  ip 

(63) 

o 
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I 


£  =  -  6  Cob  <p 


m  =  0 
z  r 


n  =  Cos  0 
z  r 


f 


J- 

1 

(64) 


(65) 


(66) 
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Figure  103  Yaw  Rotation 


Figaro  104  Roll  Rotation 


Figure  10S  Pitch  Rotation 


10. 0  STABILIZATION  SYSTEM  BRIDGE  NETWORK  DEVELOPMENT 

10. 1  Roll- Yaw  Bridge  Considerations 

The  louver  and  pitch  fan  thrust  reverser  hydraulic  actuators 
have  dual  input  valves  which  allow  the  summing  of  electrical  and  mechan¬ 
ical  inputs  at  the  same  flapper.  The  electrical  input  is  designed  so  that 
it  has  25%  of  the  authority  of  full  mechanical  input.  The  actuators  employ 
mechanical  position  feedback  to  the  valve  spool,  so  that  for  a  constant 
electrical  input  a  fixed  output  displacement  will  result.  The  valves  have 
dual  1500  ohm  coils,  and  in  the  series  connection,  4  ma  is  required  for 
full  electrical  output.  In  the  parallel  coil  connection,  8  ma  is  required 
for  full  output,  and  8  ma  will  drive  the  electrical  input  to  full  output  through 
a  single  coil.  Due  to  the  push-pull,  push-push  type  mechanism  of  develop¬ 
ing  rolling  and  yawing  moments  from  the  exit  louvers,  it  is  possible  to  con¬ 
nect  the  valve  coils  from  the  four  sets  of  louvers  into  a  bridge  circuit  such 
that  a  single  amplifier  will  drive  all  four  coils  in  the  proper  phase  for  a 
yaw  input,  and  another  single  amplifier  will  independently  drive  all  four 
coils  in  the  proper  phase  for  a  roll  input.  In  addition,  if  the  two  coils  in 
each  valve  are  connected  into  opposite  bridge  legs,  the  effects  of  an  open 
failure  of  a  single  coil  will  be  minimized. 

A  requirement  for  the  amplifiers  to  drive  the  bridge  is  that 
they  have  a  high  output  impedance.  If  they  had  infinite  output  impedance, 
i.  e.  were  current  sources,  there  would  be  no  effect  on  the  system  if  a 
single  coil  opened.  This  says  that  if  the  output  load  was  doubled  in  re¬ 
sistance,  the  voltage  would  double  and  the  current  would  remain  the  same. 

Amplifier  Input  Network  Considerations 

The  amplifier  input  network  picked  to  mechanize  the  desired 
transfer  function; 

—  —  =  K  a  .  11  *  H  T  ,  is  shown  in  Figure  106. 

P  1  +  r  S 
r  comm,  *  w 

This  approach  is  feasible  because  rate  gyros  with  high  output 
voltage  gradients  have  been  specified.  Our  rate  gyros  will  have  an  out¬ 
put  gradient  of  .  4  volts/deg/sec. 


f 


i 
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K  ^  is  varied  by  variation  of  the  amplifier  gain,  and  R  can 
be  any  of  six  values,  which  are  .  1,  .  05,  .  03,  .  02,  .  01,  and  .  0025. 

r  has  been  chosen  as  10  seconds,  nominal.  The  values  of 
t  ,  and  of  the  six  K  values  are  affected  by  the  leakage  resistance  of  the 
330  mfd.  capacitor. 

The  transfer  function  for  the  circuit  of  Figure  107,  including 
the  effect  of  leakage  is: 


Vrl 


R,  +  R0  +  Rt 
1  2  L 


RrtRT  CS 

1+^ 


R.ffl,  +  R2>cs 

1  +  R-  +  R  +  R 


L  J 


(1) 


Figure  107  Servo  Input  Network  with  Leakage 


If  we  pick  position  six  of  the  ratio  (K)  selector  switch,  R  =  3, 000  ohms, 
R t  =  27,  000  ohms,  and  C  =  330  mfd: 

thus,  writing  resistance  values  in  thousands  of  ohms, 


3  +  R- 
JL 


30  +  R_ 
L 


i  +. 


3  +Rt 

_ _ _ L_ 

30RlCS 

1+  30+Rl 


We  find  that  the  ratio  of  the  time  constants,  K,  is  then  equal  to 


3R_  30  +  R_  30  +  R_ 

Li  Li  Li 

3  +  R_  *  30  R  =  10(3  +  R  ) 

Li  L  Li 

if  we  further  let  R  =  100, 000  ohms; 

The  ratio,  K,  then  becomes 

130  _  J30_  _ 

10(130)  "  1030 


(3) 


Figure  108  is  an  amplitude  vs  frequency  plot  of  one  channel 
of  the  breadboard  stab  system  package  with  the  selector  on  position  6. 

1  18 

The  best  fit  to  this  data  indicated  that  t  =  10. 6,  and  K  =  ■  *  =  .  111. 

10.  6 


A  further  cause  of  a  ratio  error  is  the  effect  of  using  a  full- 
wave  input  demodulator  which  feeds  directly  into  the  compensation  net¬ 
work  without  filtering  of  the  A.C.  components. 

If  we  assume  a  peak  demodulator,  output  of  IV,  the  RMS 
value  of  the  output  is  .  707V.  The  DC  value  of  the  output  is  .  636V.  The 

RMS  value  of  the  AC  components  is  thus  >/.707^-.636^=  .  308V.  The 
AC  components  are  operated  upon  by  the  network  as  if  the  capacitor 
were  short-circuited,  and  are  thus  attenuated  to  a  value  of  .308K  volts. 

At  very  low  frequency,  the  DC  component  goes  through  the  network  un¬ 
attenuated.  The  RMS  value  of  the  output  at  low  freauencies  then  is 

yf,  6362  +  .  3082K2. 

At  K= .  1,  (our  max.  value)  the  RMS  output  is  equal  to 

yj,  405  +  .  00095  * .  636V.  At  high  frequencies,  both  the  AC  and  DC  values 

are  attenuated  by  the  factor  K,  (.  1  in  this  case)  and  the  RMS  value  of  the 

9  2  0707 

output  *  .  0636  +  .  0305  =  .  0707V.  Thus  the  ratio  is  =  .  111. 

instead  of  ,  10,  At  all  K  values,  this  effect  changes  the  ratio  by  the  same 
percentage.  Thus,  at  K=.  01,  the  ratio  would  be  .  0111. 

It  is  interesting  to  note  that  an  R_  of  240K  would  also  give 

Lt 

K  =  .  0111  at  the  setting  of  .  01. 
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A  Data  From  8-29-62  Freq 


Figure  108  Package  Yaw  Channel  Response,  Position  #6 


